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I Introduction 



I. 1 Curiosity and doubt 

The knowledge of mankind has grown rapidly in the last century. Especially through 
the fundamental research in physics. It was and is driven by human curiosity, since 
typically the fundamental research has no direct application in everyday or technical life. 
Nevertheless, over the time these insiglits into natural laws have shown their profitableness 
in the technical revolutions. 

So on the one hand, the growth of the knowledge of mankind has always been driven 
by curiosity, since "curiosity is always the starting point for solutions to a problem"^. On 
the other hand, physical models and theories can never be right or true in a rigourous 
sense, but only falsified. So beside our success in describing nature, "we absolutely must 
leave room for doubt or there is no progress and no learning. There is no learning without 
having to pose a question. And a question requires doubt"^. In order to learn about 
nature, we have to be curious and ask new question, but we also must have the courage 
to doubt our new found answers. 

But curiosity is not a feature of science itself, it is an innateness of mankind. Most of the 
questions asked by physicists may seem academic and irrelevant to the scientific layman. 
Maybe still complicated, the questions that are tried to be answered by nuclear phyicist, 
by colliding atomic particles at the speed of light, are those questions every human reflects: 
Where are we from and why are we here? 

Where we are from is widely believed to be the so called big bang, during which time, 
space, matter and energy turned into existence. Why or what exactly happend at time 
zero no one can tell because at this singluarity in spacetime the law of physics broke down. 
Asking about before also makes no sense, because without time there is no before. 

Why we are here is a more philosophical question which can be understood in many 
different ways. While the physicist wonders about the asymmetry between matter and 
antimatter, which is an essential fact for our existence, others would ask this question 
more spirtiually. But the creation of the universe via a big bang would not contradict 
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the existence of an omnipotent power that is responsible for "hfe, the universe and every- 
thing" i... 

I. 2 The fundamental building blocks 

But let us go back to a question we believe to have answered to a wide degree: What 
are we made of? This question has a long history and mankind has had great success in 
finding the fundamental building blocks that make up the matter around us. Most all of 
these findings where mainly made in the last century. 

The theory of the greek classical "elements" earth, water, air and fire, that go back to the 
philosopher Empedokles, persisted to the sevententh century until the new definition for 
an element by the chemicist Robert Boyle marked the beginning of the modern chemistry. 

The chemicists classified many new elements, which soon required a new ordering 
schema. The periodic system of the elements was introduced which sorts the elements 
by their mass and attributes to groups of eight. The electron, the first subatomic particle, 
was discovered in 1897 by Joseph John Thomson while its existence was predicted by 
George Johnstone Stoney 23 years before. Thomsons so-called "plumpudding"-model of 
the atom was replaced in 1911 by Ernest Rutherford after his famous scattering experi- 
ment. With Rutherfords interpretation, where a positively charged core is surrounded by 
a cloud of electrons, it was then possible to explain the chemical properties in the peri- 
odic system by the atomic electron shells. In 1919 Rutherford discovered the proton, and 
with the discovery of the neutron in 1932 by his student James Chadwick the description 
of the atomic nuclei was complete. On the one side, the atomic physics focused on the 
electron shells. Max Born postulated discrete orbits for the electrons and with the advent 
of quantum mechanics physicists were able to explain the atomic spectra. 

On the other side the nuclear physicists were studying the structure of the atomic core. 
Around the 50's more and more new particles like the proton and neutron, called hadrons, 
were found. Soon there was such an overwhelming richness, that now the nuclear physicists 
needed a new ordering schema. Murray Gell-Mann proposed his "eightfold-way", where 
the hadrons could be grouped to octets of similiar mass. This was the predecessor of 
the quark model in which the protrons, neutrons and all other hadrons are no elementary 
particles, but made up of quarks, where the mesons are quark-antiquark combinations and 
the baryons consist of three quarks. 

Today the standard model of particle physics classifies the fundamental particles into 
fermions as the particles of matter and the bosons as the particles mediating the interac- 
tions. The standard model (SM) has been very successful and marks a milestone in our 
understanding of the elementary particles and their interaction. But it falls short of being 
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a complete theory since it lacks the description of gravity. Also the many free parameters 
that have to be determined by experiment are a big shortcoming of the model. 

So there are many attempts to extend the SM (e.g. supersymmetry) or to discard it 
in favor of a new theory (e.g. stringtheory). But up to now it is the best answer to the 
question "What we are made of". 

I. 3 Back to the Big Bang 

When we go back in time, we can see different forces responsible for the evolution. At 
the beginning the universe was extremly hot and dense and all fundamental particles were 
unbound and free like in plasma. Due to the expansion the temperature was going down 
and the quarks were pulled together by the strong force and formed hadrons; this is called 
hadronization. Most of them are unstable, but protons and neutrons assembled to stable 
clusters. During the ongoing cooling the electrons were bound to these clusters by the 
electromagnetic force. The atomic elements and molecules began to develop. Finally, the 
gravity pulled the neutral atoms together and so larger clusters of matter started to grow 
forming the stars and planets making up our known universe. 

The gravity and the electromagnetic force are quite good understood, so one main topic 
of current research is the study of early universe, the regime of the strong force. In fact, 
this force is so strong that it is not possible to separate two quarks. Therefore, no one has 
ever seen (or measured) a quark. The hope is that if we could heat up the hadrons to the 
temperatures directly after the big bang, we could create this plasma state of the quarks 
where they are free. With the largest machines in the world, physicists collide nuclei at 
the speed of light to create such a hot and dense system resembling the universe a few 
nanoseconds after the big bang to study the phase transition of the strongly interacting 
matter. The theoretical description of these collisions are very involved and due to the 
complex dynamics of the fundamental degrees of freedom, this is not yet understood from 
first principles. 

In this thesis, I will analyze a phenomenologial model called recombination or coales- 
cence which is a relevant mechanism for hadronization in heavy ion collisions (HIC) at the 
highest energies that are available today. It is structured as follows: 

• In the next chapter I give a short introduction into the phenomenology HIC, explain 
different models to describe the collisions and discuss how to probe the QGP. 

• The third chapter adresses the hadronization from a QGP, focusing on the recombi- 
nation approach and its ingredients. 

• The results are presented in chapter four where I first give an overview of the pa- 
rameters and parameterizations, followed by a detailed discussion on the different 



results and predictions of recombination. 

The last chapter contains final remarks and conclusions. 
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II. 1 Phase diagram of the strong force 

The main goal of the current and past heavy ion programs is the search for a new state of 
matter called the Quark- Gluon-Plasma (QGP) [Bas99]. The quarks are the constituents 
of hadrons. They interact via their color charge by the strong force meditated by the 
gluons. Due to the strong coupling quarks can not be observed as single particles, only 
as clusters, where the color of all quarks add up to "white". This so called confinement 
is expected to break at sufficiently high temperatures and/or densities. The quarks and 
gluons are then free unbound particles like the electrons and protons in a plasma. That 
is why this state of matter is referred to as QGP and describes a new phase different to 
the normal ground state. 

The nuclear matter in the accelerators 
is highly compressed and heated up during 
the collision. At sufhenctly high beam en- 
ergies, we expect the system to go through 
a phase transition from the hadronic phase 
to the QGP. Due to the large pressure, the 
systems expands rapdily, cools down and 
crosses the phase boundary again. The 
quarks hadronize and form new hadrons. 
In this still very dense system, the hadrons 
scatter mostly inelastic and create exited 
hadronic states until the system is too di- 
lute. At this so called chemical freeze-out, 
the inelastic cross sections are neglectable. 
After the kinetic freeze-out, there are also 
no more elastic interactions between the 
hadrons. 

The livetime of the QGP is very short and only the hadrons coming from the kinetic 
freeze-out can be measured in the detectors. So there is no way to see if we have created 
a QGP or not. There are two different ways to examine the early phase of the collision: 

• We can study different probes that can escape the early system mostly undisturbed 
and try to extract information from these measures. That would be particles with 
a low cross-section in the medium like photons, pair-created leptons (dileptons) and 
heavy-quark mesons (J/ip). 
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Figure II. 1: The (assumed) phase dia- 
gram of the strong force (picture taken 
from [CBM07]). 
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• The other possibihty is to look for characterisitic properties of the hadrons from 
freeze-out, which can be associated with a QGP formation, e.g. collective flow 
phenomena like the elliptic flow. 

II. 2 Theoretical descriptions of heavy ion collision 
II. 2.1 Statistical thermal model 

The thermal model is a statistical approach to particle multiplicity in heavy ion collisions. 
One assumes a globally equilibriated thermal source described by temperature T and 

chemical potentials /7,q for the charges Q. The thermodynamical observables are then 
evaluated as an average over the statistical ensembles [BM03]. Therefore, one has to 
calculate the grand canonical partition function 



where the Hamiltonian H depends on the equation of state (EoS), which connects the 
energy density with the pressure. The EoS depends on the degrees of freedom and so 

changes throughout the evohition of the system. The EoS of the QGP is subject to 
current research and can only be determined within QCD. Below the phase boundary the 
system can be described by a hadron resonance gas with a known EoS. 

The applicability of the thermal model is very limited. The assumptions of global equi- 
librium and a static source are at best questionable. The predicted particle multiplicity, 
when using the EoS of the hadron gas, are in good agreement, but since this does not 
require the formation of a QGP, the predicitve powers are very limited. 

II. 2.2 Quantumchromodynamics (QCD) 

The only possibility to really calculate the shape and order of the phase transition would 
be to solve the equations of motion for QCD, the gauge theory for the strongly interacting 
colored particles. The approach from Quantumelectrodynamics (QED), where pertubation 
theory is very successful because of the small coupling constant a « 1/137, is barely 
applicable in QCD due to the large coupling of the strong force and the non-abelian nature 
of the corresponding symmetry group SU{?>). The pertubative QCD (pQCD) works only 
at large momentum transfers, where the running coupling constant becomes small. 

The other way is to discretize the langrangian and solve it numerically on a lattice 
(LQCD). Since this requires vast computational powers, there is no chance of a dynamical 
description from LQCD. But it can give insights into the order of the phase transition, the 
critical energy density and the equation of state [Fod02, dF06, Ste06]. These are needed 
as input to hydro- and thermodynamical models. 

Due to the lack of lattice data at high baryo-chemical potentials, I will later use the 
simple phenomenological MIT bag model to describe the confinement and the phase tran- 
sition. 



(II.l) 



11.2 Theoretical descriptions of heavy ion collision 



7 



II. 2.2.1 Fragmentation (pQCD) 

The hadron production at sufficiently large momentum transfer can be described by pQCD. 
The quark-antiquark pairs created in hard scatterings are connected by a color string and 
as they depart from each other back to back, the potential energy between them rises until 
the string breaks and creates a new quark-antiquark pair. This process continues until 
most of the kinetic energy of the original parton is converted. The hadrons, formed from 
the created quarks, carry a momentum fraction of Pa and are all emitted in the direction 
of the original parton. This is called a jet. 

The probability that the jet of parton q emits a hadron h with a momentum fraction 
z is called fragmentation function Dq^hi^) and can be measured in e'^e~ annihilations. 
The invariant cross section for a hadron h with momentum P can then calculated with 



For proton-proton collisions the situation is comparable to the e+e~ annihilations, but in 
heavy ion collisions (A-A collision) with a large fireball volume, the jets lose energy as 
they travel through the medium and will be modified. This modification is prominently 
observed in the so-called away-side jet supression: When two partons are created near the 
edge of the fireball, one jet can be emitted directly into vacuum while the other one has 
to travel through the whole dense medium and will therefore be suppressed. 

But for a dense medium, the description in terms of a jet from a single quark is ques- 
tionable and one has to account for multiple parton fragmentation (higher twist). In the 
extreme case of a very dense phase space with abundant quarks, they might simply re- 
combine. That means e.g. if a u and d quark are near in phase space, they are confined to 
form a vr^. This recombination mechanism plays an important role in the mid-pT range 
while it is dominated by fragmentation at high-pT- I will dicuss the recombination (also 
called coalescence) approach in the next chapter. 

II. 2.2.2 MIT bag model 

In the MIT bag model [Cho74, Cho75, Myh84, TheSO] the quarks and gluons are confined 
in volume with radius R equal to the radius of a nuclei. Outside of this volume, the virtuell 
quark-antiquark-pairs excert a pressure B on this bag, which prevents the partons inside 
to escape. They can be described by thermodynamics as ideal gas. 

The pressure of this gas depends on the temperature and the density which is expressed 
by the chemical potential /n^. If this pressure becomes greater than the bag pressure, the 
partons inside can escape their confinement. 




(11.2) 
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So the phase boundary ^ within the bag model is at 



1 

C 



(11.3) 



TT 



with ^ ~ i'^gq + '^gg) critical temperature Tq = 175 MeV at zero chemical 

potential. 



II. 2.3 Transport theory 

The transport model [Bas98, Ble99] is a microscopic approach which models the trajectory 
of every single particle. The evolution of the system is calulated with the Boltzmann 
equation 



df 



(11.4) 



where the collision term depends on the cross-sections and determines the time-dependence 
of the distribution function. Since only two-body collisions are considered in the gain and 
loss terms, one has to assume a dilute gas of particle with large mean free path so that 
three- or more-body collisions are neglectable. But the advantage is that since it is a 
microscopic theory, it can describe non-equilibrium systems. 

II. 2.4 Hydrodynamics 

The complement to transport theory is hydrodynamics, where one assumes a local equi- 
librium within a dense system, e.g. a very small mean free path of the particles like in 
a fluid. It is a very common model in HIC [Kol03, KolOO], especially, since the matter 
created at the Relativistic Heavy Ion Collider (RHIC) seems to behave like a perfect fluid 
(minimal viscosity) [Rom07, Tan06]. 

With the equilibrium assumption, we can set the collision term in eq. (II. 4) to zero, 
and with the baryon number current A'^'' and the energy-momentum tensor T^'', the basic 
formulas read 



d^.N'' =0, 
df.T"" =0. 



(11.5) 
(11.6) 



State-of-the-art calculations are 3-dimensional with three distinct fluid, e.g. two for the 
colliding nuclei and a third for the evolving collision zone. For simpler systems the equa- 
tions can be solved analytically which has been done by Bjorken [Bjo83] and Landau 
[Lan53]. 



1 The derivations of this formula can be found in the appendix A. 1. 
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The hydro results depend on the needed input, namely the initial conditions (Glauber 
model or CGC), the equation of state (LQCD) and the freeze-out prescription. 

II. 2.4.1 Collision geometry 

The interpretation of the time evolution of the system in the transverse plane is motivated 
by hydrodynamical ideas. The high compression of the nuclear matter generates a high 
pressure in the fireball. The difference of the inner medium and the surrounding vaccum 
results in a pressure gradient. This leads to an expansion in the transverse direction. The 
initial transverse area, determined by the overlap zone of the nuclei, increases rapidly until 
the freeze-out. So the spectrum of the particles can be described by the emission from 
a boosted thermal source, where the transverse velocity of the particles comes from the 
thermal motion with a superimposed transverse boost. 

Based on these ideas, the particle spectra are fitted by the so-called blast-wave model 
that describes a thermal source with a transverse and a longitudinal boost. Because the 
longitudinal rapidity is generally taken to be equal to the space-time rapidity by invoking 
Bjorkens longitudinal boost invariance, the model is characterized by three parameters 
namely temperature T, baryo-chemical potential fj,B and transverse expansion rapidity 
riT- 

For collisions with non-zero impact parameter the radial size and therefore the pressure 
gradient varies with the azimuthal angle (p. It is defined with respect to the plane in which 
the nuclei collide. The direction with ip = oi ip = tt is called "in-plane" and perpendicular 
to that "out-of-plane" (Fig. 11.2(a)). The protons and neutrons in the collision area that 
will take part in the scattering are called "participants", while the other nucleons, which 
just pass by, are called "spectators". 

The radial size in-plane is smaller than out-of-plane which leads to a larger in-plane 
pressure gradient. The spatial azimuthal asymmetry therefore generates a momentum 
asymmetry that is perpendicular to the spatial one (Fig. 11.2(b)). 

During the expansion, the increasing momentum asymmetry weakens the spatial asym- 
metry and so destroys its own origin; it is self- quenching. So the spatial eccentricity at 
freeze-out will differ from the initial one (Fig. 11.2(c)). Depending on the hydrodynamical 
evolution and the freeze-out time it can be smaller, zero or even negative. A negative 
freeze-out eccentricity would mean a change of direction from an elongation out-of-plane 
to in-plane. 
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■< »-l 

Impactparameter 

(a) 




Figure II. 2: Collision geometry in the transverse plane, a) Impactparameter. b) The spa- 
tial asymmetry is converted into a perpendicular momentum asymmetry due to the differ- 
ent pressure gradients, c) The initial spatial asymmetry weakens during the expansion, so 
the eccentricity at freeze-out is smaller than the initial one. 
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In the previous chapter, I shortly discussed different theories for the description of heavy 
ion colHsions and the phase diagram of the strong force. Currently, the collision dynamics 
are studied extensively within microscopical transport theory and hydrodynamics while 
the location and nature of the phase transition is subject to research in lattice QCD. 

The main question of my thesis is, assuming we have crossed the phase border in a HIC 
and created a QGP, how do the quarks form hadrons again when the plasma cools down. 
Therefore I will analyze an analytic model for the colinear recombination approach [Fri03]. 
Beside studying probes from the early phase of the collision, this offers the possibility to 
study the influence of a QGP creation on different observables at the freeze-out. 

This chapter has the following structure: 

• The first section contains the recombination formalism, 

• the second one addresses the form of the hypersurface where the phase transition 
and freeze-out take place, 

• the third one introduces the blast-wave model which models the quark densities 
depending on the collision geometry and 

• the fourth section will contain the formulas for the observables from the recombina- 
tion approach with the quark density and freeze-out hypersurface as dicussed. 

III. 1 Colinear recombination of quarks 
III. 1.1 Introduction 

The results from the Relativisitic Heavy Ion Collider (RHIC) are a confirmative sign 
for a very dense phase space so that recombination seems to play a major role in the 
hadronization process. While the strong nuclear suppression (ratio of A-A to p-p collisions) 
of pions is widely seen as the experimental confirmation of jet quenching predicted by 
pQCD, there are several key observations that cannot be explained by fragmentation. 
These are the very different results for mesons and baryons namely 

• the nuclear suppression factor, 

• the elliptic flow and 

• the large baryon / meson ratios. 
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For collisions with a higher center-of-mass energy, the influence of recombination can be 
expected to be even higher. The predictions for the Large Hadron Collider (LHC) are 
therefore of great relevance. 

I will now discuss the implications of an analytical coalescence formalism that con- 
strains the recombining quarks to have colinear momenta. For a detailed discussions of 
the derivations I refer to [Pri03]. While the onset of fragmentation sets the upper Pt 
bound of applicability, the colinearity constraint leads to an energy violation at low Pp- 
Therefore, the results for the transverse momentum spectra of this recombination model 
will be valid only for Pt ^ I GeV. But the most prominent and promising observable of 
recombination, namely the elliptic flow V2, can be expected to give valid results down to 
several hundred MeV [Pri03], since the azimuthal angle dependence in these flow coeffi- 
cients is only expressed relatively to the absolute yield. And also the very good description 
of the data (see sec. IV. 4.4) justifies the use down to low px- So such an analytical model 
can still serve as a good guidance in studying the general features and point the way for 
dynamical studies. 



III. 1.2 Non-relativistic model 

Let us start with a volume V with a homogen distribution Wa{p) of various quarks a 
which will therefore only depend on the momentum. We assume that only constituent 
quarks with colinear momentum pa will recombine to form a hadron with momentum 
Pfi = ^aPo'- This simple picture then leads us to a first non-relativistic equation for the 
hadron multiplicity which is the momentum integral over the product of the quark densities 
times the probability for these quarks to form a hadron (namely the wavefunctions overlap 
squared). For simplicity, I will stick to mesons (quark-antiquark pairs) and will give the 
generalization to baryons later. The spatial wavefunctions for the quarks and the meson 
M are: 

{x\qi,Pi;q2,P2) = y-ie^(fi-i+f2-2) (IILl) 



{x\M,Pm) = y-i/'e^(^-^^VM(xi - X2) {111.2) 

Defining the center of mass and relative coordinates resp. momentum as R= {xi + X2)/2 
and y = xi — X2 resp. P = (pi +p2)/2 and q = pi — P2, the overlap is 

{qi,pi\q2,P2\M,PM) 

=1^-3/2 j ^3 j3 ^^g^p ^- (p^ . ^ _ p^^^ _ p2a;2)] <^M(y) 

=y-3/2 y d3 i?d=^ yexp [i {{Pm-P)R- qy)] <i>M{y) 
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4^m{<i) is then the fourier transform in the relative momentum. Due to the absolute value 
squared, the deltafunction has to be rewritten as 



5\Pm -P) = j d3 xexp [i (Pm -P)x] = 

Using this in the proposed formula for the hadron multiplicity 

d^P d^pi d^p2 



Nm = CmV 



(27r)3 (27r)3 (27r) 



r Wq^ (jPl)Wq^ {jp2 ) | (gi ,pi ; ^2 , 1 M,P) | ' 



(III.4) 



(III.5) 



with a spin degeneracy factor Cm (mesons) or Cb (baryons) respectively, one finds the 
momentum distribution 



dNM 
d^P 



V f d^q P P 



Assuming an exponential parton distribution Wa{p) = To^xp 
drop from the density product. Due to the normalization of t 
vanishes and one is left with 



P_ 
T 



(III.6) 



, the q'-dependence will 



dNM ^ V 

Cm ,„_xq 7a7bexp 



d^P 



(27r) 



P 
T 



le wavefunction, the integral 



(III.7) 



The spectrum does not depend on the form of the wavefunction, it is just exponential. 
The only dependence on the specific hadron is in the degeneracy factor and the quark 
fugacities 7. 

Even in this simple approach one can explain the large baryon/meson ratio at RHIC, 
which pQCD fails to predict. It simple comes down to the degeneracies 



dNs _ 
dNM ~ Cm 



(III.8) 



which then results a proton-to-pion ratio of about 2, in contrast to about 0.2 from pQCD 
fragmentation. 



III. 1.3 Relativistic 

To make the approach more general and consistent, one starts from a density matrix Pab 
for the quarks a and b. The number of mesons is then given by 



^M = Y.i ^{M,P\p^\M,P) 



a.b 



(III.9) 
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Inserting complete sets of coordinates 

^3 



Nm = Y,J ^d3fid3f2d3fid3r^(M;P|fi,f2)(fi,f2|/9„6|fi,f^ (III.IO) 



a,b 

and using the definition of the wigner function 

{'ri,r2\Pab\r'i,f2) = J (P^'^'i + Wab{ri,r2;pi,P2)- (HI. 11) 

With the center of mass (cm.) coordinate ri^2 = (^1,2 + ^1,2) and the relative coordi- 
nate r'12 = fi,2 — r'l 2 of the correlation coordinates, one can define the relevant quantities 
for the hadron. By assuming small spatial variations in the wigner function, after some 
steps one arrives at 

^d^-^^y^ (2^0^ J (2^ 

X Wa{R; P/2 - q)^M{q)wb{R; P/2 + q) (III.12) 

with the relative momentum q between the quarks, the spatially integrated wigner function 
of the meson ^M^q) and the future orientated unit vector u'^{R) on the freeze-out hyper- 
surface S. When going to local light cone (LLC) coordinates, one can restrict the quark 
momentum to colinearity by writing it as a fraction Xa = -f^ of the hadron momentum. 
By introducing the ansatz for the light cone wavefunction 

LLC 



^(g)^irv(x) = ivn^a 



X for quarks 

\/30a::iX2 for mesons , (IIL13) 

12 \/35a;iX2iC3 for baryons 



the formula for the invariant yield becomes 

^I^^^X^W^ / D^l<^(*)l'xn^a(^;^a^) (nLi4) 

with X = {xi,X2, ..,Xnq), the number of constituent quarks ng and 

jBx = s(^l-Y,^a^Yl^^a (111-15) 



In the case of quarks this is simplified to 



E-^ =g / .3 w{R;p) (111.16) 



III. 1 Colinear recombination of quarks 
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As already shown in the previous section, the recombination formalism alone is capable 
of explaining the large baryon/meson ratio, but the results will generally depend strongly 
on the parameterizations of freeze-out hypersurface S and the used quark density Wa{R;p). 
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III. 2 Freeze-out 

Now I will model the freeze-out hypersurface. The spatial four-vector 

Xu = {t,x,y,z) (III. 17) 

will be translated into the new coordinates 

xl = {T,p,cp,v) (III. 18) 

with the transverse radial variable p, its angle (p, the time r = Vf^ — z"^ and the space-time 
rapidity 

,= iln^. (III.19) 

For the study of collisions with a nonzero impact parameter, I will directly generalize to 
an ellipsoidal cylinder. The elliptic transverse freeze-out area 



has the numerical eccentricity 



_ jR'y -Rl _ 

^numerical — ^/ (III. 21) 

if Ry > Rx- To avoid ambiguities with other eccentricity definitions, the numerical ec- 
centricity of an ellipse with the above definition will be called e while the / denotes the 
freeze-out value. 

So the surface radius depends on the angle p, which is measured with respect to the 
X-axis: 

r{p) = pf{p) with (III.22) 

N 



^Jl-ej sin2((^) 



^/a'^ sm^{p) + a"^ cos'^{p) 

\J (q!^2 -|- o?') -|- cos(2</?)(a~2 — a^) 



(III.23) 
(III.24) 



with a= (III.25) 
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The radial coordinate p is constant for a given ellipsoid surface. The length of the short 
axis X (in-plane) and respectively the long axis y (out-of-plane) are 

= r(0) = Np (III.26) 

The constant N is chosen such that the transverse is independent of the eccen- 

tricity and only depends on the parameter p. Prom 

2 

A{. = ttRA = iV^^ (III.28) 
if follows 

N = a, (III.29) 

Ra, = ap and Ry = ^ (III. 30) 

Looking at eq. (III. 23), one directly sees that the axes are changed under the transforma- 
tion a — ^ a"^ and the ellipsoid is flipped by 90°. This will be used in the next paragraph. 

III. 2.1 Integration measure 

To obtain a 3-dimensional freeze-out hypersurface, one needs to apply a constraint to the 
four- vector Xi^ or respectively x'^. This is done by constraining the freeze-out time r to be 
a function of the other three variables, so generally r = T{p,(f,r)). 

The simplest choice would be to use a constant freeze-out time with r = tq, but to be 
a bit more general and to connect the transverse with the longitudinal dynamics, I will 
perfom the derivations with an explicit p dependence only. 

According to the Cooper-Frye prescription [Coo74] the integration measure is given by 

p^da'" (III.31) 

where da^ is the normalvector on the hypersurface in the new coordinates which is cal- 
culated by a generalized cross-product of the tangent vectors. Therefore, I will calculate 
the Jacobian 

J.i = ^^ ("1-32) 
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With the dependencies 



t = t{p) cosh(ry), 
X = Pfiv) cos(99), 
y = Pfi'P) sm{ip) and 
z = t{p) sinh(?7) 

the tangent vectors read 

•='^m1 = cosh ?7,/ cos (yC,/ sin </?,|^sinhry 

Vix '-JIX2 = {0,pi-f smip + f cos(p), p{f COS ip + f sin ip),0) 
Vfi '-JixZ = (rsinhr7,0,0,rcosh77) 



where 



a ^ — 



sm{2(p). 



(III.33a) 
(III.33b) 
(III.33c) 
(111.33d) 



(III.34a) 

(III.34b) 
(III.34c) 



(III.35) 



With the fully anti-symmetric rank-4 tensor SnnXa, the infinitesimal normalvector on 
the hypersurface is calculated as 

dcr^, = e^^xa P'^dp ip^ dtp 7]'^ dr] 

= Tp (/^ cosh r],drT (^f cos ip + f' sin y?) , drT (/ sin (p — f' cos p) , sinh 77) dp dp 

(III.36) 



Similiar to eq. (III. 33), the particle momentum is parameterized as 

Pn = {niT cosh y,pT cos{(j)), pt sm{(f)), rriT sinh y) 
which then leads to the integration measure 

Pi^da'^ =dpdpdr]T{p)p 

dr 

f{(p)mT cosh(y - r/) - —pr (/ cos{p - (j)) + f sin((^ - (f))) 



(III.37) 



(III.38) 



If I had assumed a constant r independent of p, the second term would vanish and the 
integration measure would simplify to 

Piida'^ = dpdpdr]Topf^{(p)mTCOsh{y — rj), (III. 39) 

By additionally setting = 0, one recovers the measure Topdp dipdr) rriT cosh{y — 77) 



III. 2 Freeze-out 
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Model 1 Model 2 

Figure III.l: The two extreme cases for the expansion of fireball. Left: The system ex- 
pands perpendicular to the surface. Right: The system expands radially outwards. 

for a cylindrical freeze-out, since 

ef = ^ a = l ^ /((^) = 1 ^ r{ip) = p = R^ = Ry (III.40) 

III. 2.2 Surface flow 

For an ellipsoid the radial direction ip is generally not perpendicular to the surface. It is 
related to the orthogonal direction j3 by 

tan/3 = (l — Cj) tan Lp = tan p. (III. 41) 

So for an elliptical freeze-out area I consider the two extreme cases: 

• Model 1: the system expands perpendicular to the surface 

• Model 2: the system expands radially outwards independent of the surface orienta- 
tion 

Fig. III.l depicts the possibilities. When I drop the p dependence of r (eq. (III. 39)) both 
models can be related via the integration measure [Tom05] as can be shown by rewritting 
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the X and y coordinate in eq. (III. 33) using the angle /3: 



x" =pf{ip) cos ip 



cos If 



sign (cos if) 



-y/ tan^ {(p) + a~'^ 



^Va-6tan2(^) + a-2 



sign (cos /3) 




i/q;~2 sin2(/3) + cos2(/3) 
=P(/a-^a-i(/5))a^ cos/3 



(III.42) 



and similar 



=P [/a^a-i ""^ sin j3 



(III.43) 



The integration measure becomes 



dpdpdrjTop [fl^„-i{P)] mTCOsh(y-r?), 



(III.44) 



which is similiar to the one from eq. (III. 39). When renaming P ^ ip, the only difference 
is the transformed a. As already explained above, this transformation is a rotation of the 
ellipse by 90°. So the integration over an ellipse elongated out-of-plane with the system 
expanding perpendicular to the surface is equal to the integration of an in-plane ellipse 
with a radial expansion. But this exact correspondence is only true if dpT = 0. The 
consequences of this relation are discussed in section IV. 2.3.1. 



III. 3 Blast Wave Model 



21 



III. 3 Blast Wave Model 

For the parameterisation of the partonic density I take inspiration by hydrodynamical 
ideas [Sch93, Flo04] and therefore assume a focally thermalised fireball. With the four- 
momentum of the quark a and the flow four- velocity of the system u^,, the thermal 
density for fermions with the degeneracy factor C reads 

WaiR) = C 



exp [{p^u'^ - /x) /T] - 1 
exp [- {pi^u" - n) /T] 



1 - exp [- ip^u'^ - fi) /T] 

oo 

= J](-l)*^+iexp l-k {p,u^ - /x) /T] (III.45) 

k=l 

To simplify the derivations, I only use the lowest-order term which is just Maxwell- 
Boltzmann statistic. This is sufficient, since the higher order terms are neglectable as 
I have verified. 

III. 3.1 Azimuthal momentum asymmetry 

In hydrodynamics the transverse flow is generated by a pressure gradient 

which will depend on the angle for non-central collisions. The spatial asymmetry (fig. 11.2(b)) 
has the initial numerical eccentricity 



= ^M^SZES! („,.47, 



w{h) = Ra-^ and (III.48) 



h{h) 

h{h) = (I"-49) 

of an ellipse with the long axis h{b) (y-direction) and the short axis w{b) (x-direction) which 
describes the collision in the transvese plane and only depends on the impact parameter 
b. So in the x-direction there is a stronger pressure gradient as in the y-direction. This 
yields a larger expansion velocity of the system and this momentum asymmetry will be 
perpendicular to the spatial asymmetry. This expansion leads to a decrease in the spatial 
asymmetry and therefore a decelerated growth of the momentum asymmetry which is 
called a "self- quenched" behaviour, since the momentum asymmetry destroys its own 
origin. 
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Figure III. 2: Comparison of the impact parameter dependent eccentricity e from Glauber 
model and CGC [Hir06] calculations to an eccentricity from constant density (eq. III. 52) 

The size of the initial transverse area Aij, is determined by the impact parameter: 



It is not a piori clear how the spatial eccentricity will evolve until the freeze-out, but 
it generally depends on the initial spatial eccentricity, so ej = ej(ej). Depending on the 
dynamics, the system could 

• freeze-out early and leave a remaining spatial freeze-out eccentricity ej < Cj, where 
it is still elongated out-of-plane, 

• freeze-out so that the initial spatial eccentricity has been compensated by the ex- 
pansion and it is now circular with ej = 

• or freeze-out late so that the initial spatial eccentricity is overcompensated and 
reversed so it is elongated in-plane. 

The second point is just a special case of first one. So one has to distinguish between the 
elongation of the spatial area at freeze-out. In section IV. 2.3.3 I will give evidence for the 
first option and also relate Cj to e/ to study the effects of a non-circular freeze-out area. 

III. 3.2 Comparing asymmetry parameters 

I would like to state again that the used eccentricity e is the numerical eccentricity of 
an ellipse while the eccentricity definition generally used in the heavy-ion community is 



= 7rw{b)h{b) 




iTw{b)'^ 



7r{RA-b/2f 



(III.50) 
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[JacOO] 

where the average is over the transverse plane weigthed by the nucleon density. For a 
simple constant density this is simplified to 



Ry — b 

The definition of this (lets call it "geometric") eccentricity e differs from the definition 
of the numerical eccentricity e. To compare one with the other, the above equation is 
written in terms of the numerical eccentricity. 



2Ra 2 - e2 



(III.53) 



Realistic calculations of the eccentricity are the ones from the glauber model [Gla70] 
which differs from the simple approximation b/ {2Ra), because it accounts for the thickness 
of the nuclei. Especially for very peripheral collisions, where only the sparsely populated 
parts of the nuclei collide, the glauber eccentricity is much smaller. Figure III. 2 shows a 
comparison for different conditions [Hir06], including the Color Glas Condensate (CGC) 
results. The value of 6/(2i?^) is generally to large but Econst/^ = b/{4RA) is in good 
agreement with the glauber model. Therefore, I define the eccentricity by the relation 



1^ 

''22-e^ ^ ARa 



£ ■=-- 2 ^ A TD with the inverted relation (III. 54) 



(III.55) 

The relevance of this redefined eccentricity becomes visible in the following paragraph. 
III. 3.2.1 Expansion velocity 

Since the transverse flow is generated by a pressure gradient and I will not model the 
pressure profile, I take the ansatz from eq. (III. 46) and assume rjT ~ with the mean 

fm 

transverse flow rapidity free parameter. To have a much more handy expression. 
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I rewrite the denominator and expand it to fourth order in e: 
1 _ \/l-efsin2((^) 



4e 

0.5(1 -cos(299)) 



2e + l 



1- 



2e + l 



=:l + fjT{v) (III.56) 

So the elhptic asymmetry cos(2(^) scales mainly linear with the defined eccentricity 

e=^™i, (III.57) 

which is comparable to the glauber eccentricity from central to mid-central collisions as 
can be seen in Figure III. 2. This leads to a direct linear relation between the elliptic flow 
and the eccentricity which is also expect from hydrodynamics (see section IV. 4.2) and 
therefore justifies the choice in eq. III. 54. 

Together with a radial profil in p, the expression reads 

titM = VrNd ' (1 + fjri^P)) ("1-58) 

dp + 2 

The normalization is chosen to be -^-^ — so that 

Now the four velocity of the expanding fireball is created with a longitudinal boost 
rjL and a transversal boost ijt- Taking Bjorkens boost-invariance argument, I set the 
longitudinal rapidity equal to the space-time rapidity, i.e. r]L = rj: 

Uu = (cosh cosh ryr, sinh ryr cos sinh r]T sin /?', sinh r]L cosh t^t) (III. 60) 

The angle (3' depends on the model chosen for the surface expansion from section III. 2.2 
with (5' = P{<f) for model 1 (perpendicular expansion) and P' = (p for model 2 (radial 
expansion) . 



III. 4 Observables from recombination 
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III. 4 Observables from recombination 
III. 4.1 Invariant yield 

Now I combine the recombination formalism with the discussed freeze-out hypersurface 
(sec. III. 2) and the quark density from the blast- wave model (sec. III. 3). With eq. (III. 16) 
the invariant yield for the quarks read 



9 



d'^p prdpTdcpdy 

rpo p2-K 

dp j dip I dr]T{p)p 



9 



(2vr)3 



dr 



f {ip)mT cosh(j/ - r/) - —px (/ cos(/?' - (j)) + f sin(/3' - 0)) 
X exp [- {p^u^ - n) /T] 
where the exponent expands to 

Puu'^ = mx cosh(y — rj) cosh r]T — Pt cos(/?' — (p) sinh ijt 



(III.61) 



(III.62) 



To study the dependence on the particles momentum angle (f), one uses the fourier expan- 
sion 



dN 



Pt dpT d<pdy 



vo + 2 cos{n(p)vn + sin(n^)'U„ 



n=l 

dN 



2'KpT dpTdy 



1 + 2 Vn cos(n^) + Un sin(n^ 



n=l 



with 



= ^ J d(j)COs{n(j)) 



dN 



1 



Pt dpT d(f)dy 
dN 



'^n = TT I d^Asin(n0) — 

ZTT ./ pTdpTdcpdy 



The ^-integrated yield vq 



dN 



(III.63) 

(III.64) 
(III.65) 

describes the overall transverse momentum 



27TPTdPTdy 

dependence, while the other fourier coefficients Vn and tt„ encode the deviations from an 
azimuthal symmetric yield. To have a better comparison, the normalized fourier coeffi- 
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cients Vn and Un, also called flow coefHcients, are introduced: 



Vn '■=-?- = (cos(n(^)) and 
Un ■.= ^ = (sin(ra^)) 

Uq 



(III.66) 
(III.67) 



The full expression reads ^ 

KiInP{v)mT cos(n/3') 



- k J"^ ^^^"-^^ f{ip)pTdpT cos(n/3') 



K f {v)PTdpT — 



c? cos((n - 2)/?') - cos((n + 2)^') 



(III.68) 



For the coefficients ■u^, one has to replace the cosine terms with sinus. From that the 
0-integrated yield follows as 



dA^ 



=5exp [/x/T] ^l^^l^ p'r{p')p' 
KilQp{(p)mT - KoIif{ip)pTdpT 



(III.69) 



The arguments for the modified bessel functions = {k{p,ip,pT)) and Im = Im {iip^f^PT)) 
are 



j^jQuark (^p^ip^p^'^ — — T. Vri^jPjPT) respectively 



r-i quark/ x siuh //^(y^.p.pr) 

br {p,^,PT) = 



(111.70a) 
(111.70b) 



Due to the symmetry of the system, most coefficients vanish. The only ip dependence is 
in the function /(</?) which only contains terms proportional to cos(2n(^), as can be seen 



1 How to solve the integrals can be found in the appendix A. 2. 
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in eq. (III. 56). Therefore, the expectation values (sin(n0)) and (cos((2n+ 1)^)) vanish: 

U^ = = Un (III.71) 

V2n+1 = = U2„+l (III.72) 

This results in the symmetries 

dN dN 

-W= . and (III.73) 



prdpTd^dy p-rdpTd^dy 

dA^ dN 
, . , , W = , , , , - 0) (111.74) 
Pt dpT dcpdy px dpx dcpdy 

For the hadrons the integration over the product of the quark densities involves an 
additional integral over the momentum fractions x = (xi,..,x„^) times the wavefunction 
\iP{x)\' (eq. (III.14)): 

dN f Puda 



Pt d Pt d(pdy 



Due to the product Y\g^Wa{R;XaP), the arguments of the besselfunctions are now x- 
weigthed sum of the above arguments (eq. (III. 70)) for the quarks: 



[k] 



hadron/ sr^ V'ml + {xnPTf coshr]T(,V,p,XnPT) mi^ya \ 

{p,(p,x) = 2^ , (III. 76a) 



n=l 



ihadron. -.x \^ XnPTsmhr]T{(p,p,XnPT) mT7au\ 
\ {p,'P,x)=2^ Tf; (III.76b) 



71=1 



where Uq is the number of valence quarks. 

The degeneracy factor C are not a priori clear from QCD. In principle every quark has 
3 color and 2 spin degrees of freedom. But since there are no dynamical gluons in this 

model it would be consistent to require that the quarks have the right quantumnumbers. 
Therefore, the degeneracy is only determined by the degrees of freedom of the hadron, e.g. 
Cp = 2 and = 1. I will not take into account feeddown from resonance decay, except 
for the A, where the is too close in mass to be suppressed. Hence I take Ca = 4 [FriOS]. 
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III. 4.1.1 Central collisions 

For collisions with zero impact parameter, the initial as well as the freeze-out eccentricity 
will be zero. Then /(</?) = 1 and tjt is independent of (p and so one arrives at 

C^^^[p/T]^ j\p'T{p')p' j Dx 



2nPTdPTdy 



X 



KiIouiT - KohpTdpT 



(III.77) 



III. 4.1.2 Peripheral collisions 

The parameter po measures the size of the transverse freeze-out area Aj, = irp^ (see 
eq. (III.28)). 

The total multiplicity depends on the initially created quark density, which I assume to 

scale linear with the size of the initial transverse area Aj,{b) = '""^(^) (III. 50)), 

Vl - ei(6)2 

which only depends on the impact parameter b. In order to have only one parameter for 
all centralities, I let po scale like 

Po{b) = Pc ■ (III.78) 

and will fix pc for central collisions. 

III. 4.2 Flow components 
III. 4.2.1 Elliptic flow 

The first non-vanishing fourier coefficient in the 4> expansion is the so called elliptic flow 
V2 which will therefore dominate the expansion. It measures the momentum asymmetry 
between the particles emitted in-plane (x-direction) and out-of-plane (y-direction): 

V2 := (cos(2^)) = (cos^((/)) - sm^{(j))) 

For dpT = 0, it can be calculated by 

. X fn^dy? LMpV'Ki 72/2(99) cos(2/3') 

V2(Pt) = 1 —-^ (111.80) 

Cd^J,Up'p'K,IoP{^) 

The importance of this observable comes from the fact that it is generated by the self- 
quenching momentum asymmetry (see sec. III. 3.1) and is therefore mostly sensitive to 
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the initial (partonic) stage of the colHsion. 



III. 4.2.2 Hexadecupole flow 

To study additional effects, I also take a look at the next order V4. By considering the 
ratio of the different coefficients, one can disentangle the different contribution to the 
flow. Since the coefficients scale like V2 ^ s and ~ e^, one expects the ratio of the flow 
coefficients to be independent of the initial eccentricity. In the section III. 4.3.2, I will give 
a more detailed description of this ratio. 



III. 4.2.3 High-pT flow 

The flow components calculated with recombination are monotonically rising with pT- 
Since I will not consider contributions from fragmentation which dominates over recombi- 
nation at high pt, the experimental observed drop at a finite pr has to be modeled by a 
phenomenological factor [Fri03] 

t^iPT) = ,^1^, ^2 ■ 

1 + {pt/poY 

with the additional parameter pQ. It enters the transverse fiow rapdity (eq. (III. 58)) and 
makes sure that faster partons do not feel the expansion asymmetry that much: 

riT{'P,p,PT) = VtNcI (^-^^ (l + ?7r((^)Ac(pr)) (111.82) 



III. 4.3 Constituent quark number scaling 
III. 4.3.1 Elliptic flow scaling 

One major success of recombination is a predicted simple connection between the quark 
elliptic flow and the hadron elliptic flow which gives strong evidence for a QGP phase 
in the collision. This connection is called "constituent quark number scaling" (CQNS), 
because the hadron elliptic flow scales with the number of constituent quarks. 

To study the connection of the elliptic flow between quarks and hadrons, let us fourier 



30 



III Hadronization from a QGP 



expand the quark density similar to eq. (III. 63) as 

Wq{pT,'P-i4') =— J dcpWq + — '^^cos{n4>) J d (f) cos{n(f))wq + sm(n(j)) J d (j) sm(n(f))wq 



n=l 

=exp [k{ip,p,pT) cosh(r7 - y) + ji/T] Iq 



n=l 



" cos(n/3') cos{n4i) + sin(n/?') sm{nq 



oo 



=exp [kiip,p,pT) cosh(r? - y) + p/T] 1^2 ^ v^^ cos(n(/3' - 4>)) (III.83) 



n=0 



by defining := 1/2 and < := 4/^0- 

To simplify the derivation, I take dpT = which is only a minor simplification compared 
to the following ones. Because of the orthogonality of the cosine, only the coefficient V2 
enters the elliptic flow for quarks: 



r .quark /p P' d p' C ^ ^ f {^)hK,v^ cos{2(3') 
Jo ^ dp' Jo d(^/^(<^)/oKi 



To compare it to the hadron elliptic flow, I simplify the hadron wavefunction. Since it is 
the product of the momentum fractions x^, it is already maximal for quarks with equal 
momentum. Therefore, I use the delta-function approximation 



ft^fxi-— ) (III.85) 



Then the x-integration will brake down and all quarks have the equal momentum fraction 
l/n^. For a meson one now has a product of two fourier series which can be expanded to 

Wq,{pTl2)wq^{pTl2) =exp [A;"^™(v?,p,l/2) cosh(?7 - y) + 2///r] (i((^,p,pT/2)) 

00 

4 «cos(n(/3'-,^))cos(m(/3'-<^)) (III.86) 



00 

X 

n,m=0 



As I will show in section III. 4.4, the higher terms can not be neglected, so analytically 
there is no simple and exact relation between the quark and the hadron V2. 

Those who want to derive an analytical expression for the scaling, must use a strong 
simpliflcation and drop the spatial correlations between the angle and flow velocity from 
the quark density. In this case, one assumes just a constant asymmetry in the flow proflle 
which is equal at every point in the flreball, and rewrite the density with the help of the 
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quark elliptic flow as [Fri03] 

w,{pT,V,(t>) = w,{pT,(t>) (l + 2 cos(2</))) 



(III.87) 



The azimuthal asymmetry is only contained in the cos (2^) term which depends only on 
the hadron emission angle (f). The previous dependence of rir on the spatial angle (f is now 
gone. When assuming a circular freeze-out area (e/ = =^ /(ip) = 1), the (/^-dependence 
can be integrated out completely which leads to a much simpler equation for the hadron 
elliptic flow. One finds for mesons 



[V2] 



meson 
Pt 



p' d p'Iq [i(p,PT/2)] Ki [k{p,PT/2)] 
' p' d p'lo [i{p,PT/2)] [k{p,PT/2)] 



q2 



1 + 2vrv: 



(III.88) 



and for baryons 



r ibaryon 



jl p'dp'h [Kp.PT/m Ki [k{p,PT/3)] 

p'dp'h [i{p,PT/3)] Ki [fc(/9,Pr/3)] 



„93 



(\^2{vfvf^vtvf^vfvf)) 



(III.89) 



where v\ is to be evaluated at Pt/2 for mesons and at Pt/3 for baryons. With the 
additional simplifications of a constant transverse fiow rapidity independent of the radius, 
one arrives at 



r iin 



meson 



V.- 



q\ 



//2 



1 + 2vfvf 



(III.90) 



and 



Nbaryon 
Pt 



Pt/2 



ql g2 g3 



1 + 2 (vtvf 



q2 q3 



(III.91) 



Pt/3 



respectively. For hadrons with similiar quark content all coefficients i;| are equal and by 
further neglecting quadratic and cubic terms, one obtains the much celebrated constituent 
quark number scaling [FriOS] 



(III.92) 
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Beside all the simplifications, this scaling law is confirmed with great success by the 
experimental data. A common way to show the good agreement is to scale the data of 
hadrons by l/n^ and plot it as Pr/uq vs. vijuq. According to the scaling law, all curves 
for mesons and baryons will lie on one universal curve which would be the quark elliptic 
flow. 



III. 4.3.2 Hexadecupole flow scaling 



With the above simplifications, one can also derive a scaling law for the hadron v^. Again 
neglecting higher powers and assuming equal V2 and for all quarks, one finds 



r 1 meson 



ql q2 
2 ^2 



= [2V4 + ' 



,21 quark 

^2jpj,/2 



(III.93) 



and 



baryon ^ qi 



'2 ^2 



q2 



'2 ^2 



(III.94) 



Combining the both scaling laws, ratio between V2 and can be expected to be ap- 
proximately 



H meson 



1 1 
4 + 2 



l''2 



- quark 

Pt/2 



(III.95) 



and 



VA 



-1 baryon 



2J Pj, 



1 1 

3 + 3 



VA 



L^2 



-| quark 
Pt/3 



(III.96) 



III. 4.4 The breaking of the CQNS 

Although the CQNS is experimentally well observed, the scaling is broken by the masses 
of the hadrons and the simplified equations from the upper section ca not account for that. 
Therefore, I want to emphasize that the scaling laws will only serve as a rough guide and 
all calculations will be done with the full spatial correlations of the quark density, if not 
specified otherwise. 

Let us shortly discuss the scale breaking terms: Taking the quark-antiquark density 
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from eq. (III. 86) the meson elliptic flow can be calculated via 



[V2\ 



meson 

Pt ~ rl J r27r 



(III.97) 



Vq and V2 are the contributing terms from the infinite sum. Because of the orthogonality 
of the cosine, all terms not containing the integral weightings (i.e. cos(0(/)) for the denom- 
inator and cos(2(^) for the numerator) have vanished after the 0-integration. With the 
trigonometric product 

cos(na;) cos(ma;) = ^ (cos(n — m)x + cos(ra + m)x) (III. 98) 

an infinite sum of coefficient products v^v^ that fulfill n + m = d or \n — m\ = d, 
contributes to 




2 j d(^cos^(#). (III.99) 



(Ill.lOO) 



(III.lOl) 

In fig. III. 3 you see the comparison of the full calculation to zeroth order. In zeroth order, 
the terms contain only the V2 contributions: 

V2 = 2'k{2v^) Vo = 2'!t{1 + 2v^v1^). (III.102) 

This would correspond to the CQNS. Obviously, the simple approximation of zeroth order 
is not sufficient when calculating with the full momentum space correlations. Therefore, 
the higher terms can not be discarded and a simple, analytical form for quark number 
scaling does not exist in this framework. 

Additionally to the contribution from the higher order flow components, the bessel 
functions give further deviations from the scaling law: In the high Pt region the additional 
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> 




Pt 

Figure III. 3: Study of the CQNS: The eUiptic flow of pions from the full calculation com- 
pared to zeroth order (CQNS) in the expansion of eq. III. 99 together with the light quark 
flow. 



powers of Iq give an enhancement for the elliptic flow, so it will be higher than nqV2{PT / nq) . 
The higher Uq the higher the enhancement will be. On the other hand, in the low Pt 
region where Iq is approximately one, the scaling will be broken by Ki. The argument 
k^^^^^^{if,p,l/nq) is the sum over the transverse masses of the quarks at Px/nq. So for 
Pt — > the transverse mass is equal to the quark mass and the argument is proportional 
to the sum of the quark masses. Since Ki is a monotonic decreasing function, hadrons 
with heavy quarks have a lower elliptic flow at low Pt as compared to hadrons with light 
quarks. Similiarly, baryons will have a lower elliptic flow than mesons when they have 
quarks with comparable masses. 

This explicit mass scaling is experimentally observed and can not be explained within 
the simplifications for an analytical CQNS. 
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In this chapter I will present my results from recombination with the discussed quark 
density and freeze-out hypersurface. The chapter is organized as follows: 

• Before I will present the results, I would like to summarize the parameters in the 
model and how they are determined. 

• Then I will show some effects of the parameters on the observables. 

• The third section will contain results of the transverse momentum spectra. 

• And in the last section, I will discuss the flow coefficients V2 and 'U4. There, I study 
the ratio of these coefficients to extract reasonable parameter values and then discuss 
the pt, centrality and yfs dependence. 

IV. 1 Overview of the parameters 

IV. 1.1 Recombination parameters 
IV. 1.1.1 Quark masses 

The bare quark masses (current quark masses) for the light quarks are much smaller 
then the hadrons they constitute, with m„ = 1.5 — 3.3 MeV, = 3.5 — 6.0 MeV and 
vris = 70 — 130 MeV. So while the mass of u and d make up only about 0.15% of the 
proton mass, the main contribution comes from the gluons and the virtuell quark-antiquark 
pairs. To account for the dynamically generated mass, the quarks are assigned so called 
constituent masses. 

Since in the recombination approach there are no dynamical elements Uke gluons or 
pair production, it would be consistent to assume that the recombining quarks are already 
surrounded by this virtuell cloud. Therefore, one has to use constituent quark masses. 

These constituent masses are chosen to fit the hadron masses calculated from mass 
formulas within the quark model. The best fit is achieved when using different masses 
for mesons and baryons since they have different nuclear size. I take the values from 
Gasiorowicz [Gas81] as 



rriu =md = 310 MeV 
niu =171(1 = 363 MeV 



rus = 483 MeV 
nis = 538 MeV 



for baryons. 



for mesons 



(IV.l) 
(IV.2) 
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Mesonic systems consisting of heavy quarks like charm or bottom can be described 
non-relativisticly due to the large mass of the quarks. That is why I take the constituent 
masses as the half of the meson masses. That means 

rric =mj/^/2 = 1.548 GeV (IV.3) 
rrib =mY/2 = 4.730 GeV (IV.4) 



IV. 1.1.2 V\\g\\-pT damping 

To account for the fact that the elliptic flow has some maximum and is then slowly droping 
at higher Pt, I introduced a phenomenological damping factor in section III. 4.2.3. It will 
be chosen to fit the high-pT elliptic flow data. 

Since a non-zero freeze-out eccentricity will give an enhancement to the high-pT elliptic 
flow, the value for po will be coupled to ej. 



IV. 1.2 Freeze-out hypersurface parameters 
IV. 1.2.1 Transverse freeze-out area 

The transverse freeze-out area At = vrpg effects the total multiplicity. I flx po = H fm to 
fit the experimental data on the invariant yield of different hadrons in central collisions 
(see Fig. IV. 14). As already noted, the multiplicity in peripheral collisions is assumed to 
scale with the initial transverse area which depends on the impact parameter. 



IV. 1.2.2 Impactparameter 

The impact parameter h is not a real free parameter, but there are some uncertainties 
in relating it to the centrality which is the experimentalists impact measure. To do it 
without additional calculations from the Glauber, I take the functional form [Bro02] 

c{N) = with o-jnei = 7.05b = 705 fm^ (IV.5) 

''^inel 

where c{N) is the centrality of events with a multiplicity higher than N. Despite its 
simplicity, the results are almost equal to the Glauber model for binary collisions. Only 
for ultra peripheral collision {b > 14 fm) there are small deviations. 



IV. 1.2.3 Freeze-out eccentricity 

While the initial eccentricity is fixed by the impact parameter, the eccentricity at freeze- 
out, which originates from the initial one, depends on the expansion dynamics. Due to 
the induced momentum asymmetry, the eccentricity of the fireball will become smaller as 
it expands. But there is no need to expect it to be zero at freeze-out. 
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The numerical eccentricity Cf models the shape of the transverse freeze-out area. The 
consequences from a non-zero e/ are discussed in section IV. 2.3.1 and the functional form 
of the dependence on can be found in section IV. 2.3.3. 



IV. 1.2.4 Time dependent hypersurface 

By assuming an elliptic freeze-out in the previous section, the transverse part of the 
hypersurface will be treated quite general. And the assumed boost invariance of the 
longitudinal expansion in section III. 2 is generally accepted and well established. 

An open question concerns the time dependence of the hypersurface. As already said 
in section III. 2.1 I let r be a function of the radial coordinate p. The functional form of 
the dependence and its implications are discussed in section IV. 2.4. 

Together with po, the mean freeze-out time 

TO = (IV.6) 

will enter only as an overall normalization factor. Therefore, I take it to be constant with 
To = 5 fm [Fri03]. 



IV. 1.3 Blast wave paranneters 

The parameters for the blast wave model are the temperature T, the baryo-chemical 
potential pB, the transverse flow rapidity 7]^ and its radial profile parameter dp. The first 
two depend on the phase boundary and will be determined within the MIT bag model, the 
rapidity has to be extract from fits to experimental data and the profile will be adjusted 
to a reasonable value. An additional strange fugacity 7s = 7s = 0.8 is introduced to fit 
the invariant yields (sec. IV. 3.1) and ratios (sec IV. 3.3) of the strange hadrons. 



IV. 1.3.1 Phase boundary 

The thermodynamic quantities T and ps at the phase boundary are calculated via the 
MIT bag model from sec. II. 2.2.2. To choose reasonable values for the different colliders, 
one needs the dependence on the center-of-mass (cm) energy ^/s for either T or pe- 

One can extract the values for these quantities within the tliermodynamical model by 
fits to experimental data. Prom the values in [Cle06] I choose the ansatz 

T(Vi) =Tc (1 - (Vs/bf) (IV. 7) 

with a = -0.85, b = 1.33 GeV, c = 1.12 GeV and d = 5.48 GeV. The fit is shown in 
fig. IV. 1. 

But since these are the values at freeze-out, one has to link the baryo-chemical potential 
from there to the phase boundary. By applying an isentropic expansion within a simple 
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Figure IV. 1: Temperature T and baryo-chemical potential /is 

a) T and /is at freeze-out as a function of ^/s for different colliders. The lines are fits from 
eqs. (IV.7) and (IV.8). 

b) Comparison of the phase boundary (solid) from the MIT bag model and the freeze-out 
curve (dashed) from a fit to different colliders. The dotted lines are the paths along an isen- 
tropic expansion. 



hadrongas model, I establish a connection of the conditions from the phase boundary (pb) 
to the freeze-out (fo). 

The isentropic path of the system in the T — /i^-plane is shown in fig. IV. 1(b). It shows 
that the baryo-chemical potentials can be simply related by 

/ig^ = 0.938 • (IV.9) 

Together with eq. (IV.8) and eq. (II. 3) I can express T and /is at the phase boundary as 
a function of ^/s. 



IV. 1.3.2 Transverse expansion 

The mean transverse flow rapidity r/^, will be fitted to flow rapidities extracted from the 
experimental data [Xu02] at kinetic freeze-out. The data showed in [Xu02] is for mean 
transverse flow velocity /3t, but I use the rapidity -qT = tanh/?^ to avoid that the velocity 
becomes greater than c in the fit. 
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Ecm. [GeV] 

Figure IV. 2: Parameterisation of the transverse flow rapidity 77T (full line) as function of 
y/s. The data at kinetic freeze-out (crosses) are taken from [Xu02]. 

To fit these values I choose 

rj^^^^^iV^) =a + bx + cx^ + d\n{x) (IV. 10) 

with X =ln 

with the constants a = 0.418, b = —0.064, c = 0.012, d = 0.170. As these values are 
extracted at freeze-out, I scale the obtained transverse rapidities by a constant factor k = 
0.85 to obtain the transverse flow at the hadronization surface. Using these parameters, 
the value for vt = 0.54 at ^ = 200 GeV (RHIC) agrees with the value from [Fri03]. 

For ^/s = 5.5 TeV (LHC) I obtain a transverse flow velocity of vt = 0.75 also in line 
with previous estimates [Fri04]. Fig. (IV. 2) depicts the fit (line) and the available data 
on r]T (crosses). 
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IV. 2 Influence of the parameters and parameterisations 
IV. 2.1 Delta-shaped wavefunctions 

0.5 



0.4 
0.3 
0.2 
0.1 


0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 
Pt [GeV] 

Figure IV. 3: Relative deviations rg of delta-shaped wavefunctions from realistic light cone 
wavefunctions in eq. III. 13. 

The ansatz in eq. (III. 13) for the local light cone wavefunction with the product of 
the momentum fractions x is motivated by the asymptotic form of the pertubative pion 
distribution amplitude [Pri03]. This function is very broad and so the quark momenta are 
strongly smeared around P/uq. 

The other extreme case is the delta-function approximation of eq. (III. 85) which is 
already used in the derivation of the CQNS, where the hadron consists of Uq quarks 
with all having the same momentum P/uq. The influence of this strong simplification 
compared to the realistic light cone wavefunctions is shown in Fig. IV. 3, where I have 
given the relative deviations = [dN^ — dN)/dN for pions and protons. 

IV. 2.2 Effects of the blast wave parameters 
IV. 2.2.1 Baryo-chemical potential 

The influence of the baryo-chemical potential //_b on the results is quite neglectable. First 
of all, it only enters as a factor for each quark in the invariant yield. Second, it drops 
out of the flow coefficients completely at least for boltzmann statistics. Finally the value 
oi at RHIC and higher energies is much smaller then the temperature (Fig. IV. 1(a)). 
Therefore the fugacity 7 = exp [fis/T] ^ 1 gives only small corrections. 

IV. 2.2.2 Temperature 

The temperature T sets the slope of the thermal quark spectrum and also enters the 
invariant yields. So the temperature is an important parameter, but it varies mainly 
at low cm. energies. For ^/s > 10 GeV it quickly approaches the critical temperature 
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Tc = 175 MeV and for RHIC energies and beyond it is essentially constant (Fig. IV. 1(a)). 
So for the energy range, where recombination can be considered to play an important role, 
the energy dependence of temperature is not significant. 



IV. 2.2.3 Transverse flow rapidity 

The main parameter with the biggest impact on the results is the transverse flow rapidity 
r]T = atanh(/?T), since it is expected to show strong variations as the fit from fig. IV. 2 sug- 
gests. Although the functional form is very uncertain, the general trend of a rising rapidity 
can be expected. The exact shape of the dependence is therefore not that important, as 
long as it fits the value at LHC which is generally expected to be vt = 0.7 — 0.8c. 

The influence of tjt can be studied by looking at the slope of the pr spectra. The bessel 
functions (eq. (III. 69)) can be expanded in exponentials exp [—itt-t/Tqq] with an effective 
temperature that sets the slope. The inverse temperature can then be found by 



drriT 



■In 



dN 



Pt dp2 



For central collisions, a fixed radial coordinate p and dpT 
calculated analytically [Sch93]: 



(IV.ll) 

0, this expression can be 



1 

T 



drriT 



In 



rriT cosh rjT 



Pt sinh tjt 



rriT 



1 ^ Ii rriT sinh rjr 
rriT Iq PtT 



Kq cosh r)T 




rriT-^oo cosh r^T 



So at high rriT (or pt) the effective temperature is blue shifted as 



(IV.12) 



^eff = ^1 



1 + /?T 



> T 



(IV.13) 



so the slope at high transverse momenta is less steep for a non-zero transverse expansion 
velocity /3t- 

At low rriT the case is less simple. Therefore, I show the invariant yield for pions and 
protons in Fig. IV. 4 for different values of /3t and also for two different temperatures. The 
same comparison for the elliptic flow is depicted in Fig. IV. 5 



IV. 2.2.4 Radial profile of the transverse rapidity 

The second parameter in the transverse rapidity rjT, beside the mean rapidity r/^., is dp 
which models the radial growth in eq. (III. 58) as {p/po) In the derivation of the CQNS, 
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Pj [GeV] P-r [GeV] 

Figure IV. 4: The invariant yield of pions and protons for different transverse flow veloci- 
ties Pt- 

I already used the appproximation of a radially constant transverse rapidity i.e. dp = 0. 
This is a very unrealistic assumption since there is no pressure gradient in the center and 
so there should be no expansion. Thus, one expects r/y to be a smoothly rising function of 
the radial coordinate with a dependence somewhere between square-root and quadratic. 
Therefore, reasonable values for dp should be between 0.5 and 2. 

The yield and elliptic flow of pions and protons are compared for different values of dp 
in Figs. IV. 6 and IV. 7. For the pion yield the modifications are neglectable, but at low px 
the proton yield shows a strong enhancement for higher values of dp. On the other hand, 
the elliptic flow decreases with increasing dp for both particles at low to mid px- 

The real dependence can be quite non-trivial and can only be studied in dynamical 
calculations. To have a possibilty to determine a reasonable value, I will compare the 
results for the yield to experimental data (section IV. 3.1). The elliptic flow in fact is 
also very sensitive to the choice of dp, but it also depends very much on the freeze-out 
parameters as can be seen in the next sections. Therefore, a fit to V2 data would be too 
involved and uncertain. That is why I will use dp = 1, since it is the best fit for the yield 
and also a reasonable compromise between square-root and quadratic dependence. 
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Figure IV. 5: The elliptic flow of pions and protons for different transverse flow velocities 
/3t. 




Figure IV. 6: The invariant yield of pions and protons for different radial profiles of the 
transverse rapidity. 
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Figure IV. 7: The elliptic flow of pions and protons for different radial proflles of the 
transverse rapidity. 



IV. 2.3 Separating flow and non-flow effects 

The information about the (^-dependence of the invariant yield is stored in the fourier co- 
efficients Vn- There are two different mechanisms that generate the azimuthal asymmetric 
particle yields. So before going into a detailed analysis of the results from V2 and v^, I 
will analyse the respective strength of both contributions. The fourier coefficients Vn can 
be separated in two different contributions: 

• The flow contribution comes from the asymmetry in the expansion velocity which 
is parameterized in eq. (III. 56). The asymmetry depends on the initial numerical 
eccentricity e,. This is the contribution which results in the CQNS and therefore 
the more prominent one. 

• The geometrical contribution comes from the asymmetry of the transverse freeze- 
out area which depends on the numerical freeze-out eccentricity e/. Because of its 
spatial, geometrical origin it is nearly independent of the particle species. 

Both effects develop from the the same origin, namely the initial eccentricity in non- 
central collisions. While the strength of the flow part is fixed by the impact parameter, 
the geometrical contributions depend on the functional dependence e/ = e/(ei). 

Hence, I will first discuss the geometrical and then the flow contributions and after that 
study the freeze-out eccentricity. 
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IV. 2.3.1 Geometrical contributions 

To focus on the geometrical contributions, I turn off the flow contributions by setting 
Cj = 0. Without an initial spatial anisotropy, there is no expansion velocity asymmetry. 
This normally also means that the freeze-out eccentricity is zero, so I will set it manually 
to a finite value. 

As already explained in section III. 2.2 for an elliptic freeze-out, the radial direction 
is not equal to the one perpendicular to the surface. Thus, we have to investigated the 
contributions for both models since non should be excluded a priori. Fig. IV. 8 shows V2 
and tj/i contributions for pions and protons in both scenarios. Due to the geometrical 
nature, the contributions are nearly equal for pions and protons, so I will not show any 
calculations for other particle species. 

What can be seen directly is that the is equal for both cases, while the V2 switches 
sign when going from perpendicular (model 1) to radial expansion (model 2). This is 
the consquence of the transformation behavior of the integration measure derived in sec- 
tion III. 2.2. As stated therein, the flow of model 1 is equal to the flow from model 2 with 
the ellipse rotated by 90°. Reversing this rotation is just switching the x- and y-axis. For 
V2-, which compares the flow of x- to y-direction, this is equal to changing the sign, but 
U4, which compares x- and y-direction to the diagonal direction, is not affected by this 
rotation. 

To distinguish between the models, first one has to predict the orientation of the freeze- 
out ellipse. Then, by comparing the geometrical together with the flow contributions to 
experimental data, the correct model can be identified. On the other hand, this means 
that without the knowledge of the ellipse' orientation, there is no way to exclude one of 
the models. 

In section IV. 2.3.3 I will present experimental evidence for the direction of the elonga- 
tion. 

IV. 2.3.2 Flow contributions 

The flow contributions are studied by setting e/ = 0. These contributions are sensitive to 
the quark content and are the source of the CQNS discussed in sec. III. 4.3.1. To simplify 
the comparison to the geometrical contributions, I also set dpT = 0. The modifcations 
arising from r = t{p) will be discussed later. Fig. IV. 9 shows the scaled elliptic and 
hexadecupole flow for different hadrons. 

The general trend of the CQNS is clearly visible which is one of the major successes of 
recombination. Despite the simplifications made in section III. 4.3.1, the CQNS is at the 
heart of the formalism and therefore persists even when calculating with the full space- 
momentum correlations. Although the very simple scaling law is nice, it should not be 
taken to far, because the scaling is explicitly broken as discussed in section III. 4.4. Most 
obvious is the additional mass scaling where heavier particles have a smaller and even 
negative flow coefficient. 

The mass scaling between mesons and baryons comes from the sum of the constituent 
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Figure IV. 8: Purely geometrical contributions to the elliptic flow of pions (red curve) and 
protons (green curve) from an elliptical freeze-out area with ej — 0.62. 
Left: The system expands orthogonal to the surface, which yields a positive V2 (upper 
curve) and a positive (lower curve). 

Right: The system expands radial. The V2 (lower curve) is negative, but with same ampli- 
tude, V4 (upper curve) is equal to the other case. 

masses as explained in section III. 4.4. And the mass scaling among different mesons/baryons 
comes from the quark flow (Fig. IV. 10) where only the different masses lead to a different 
elliptic flow. This is also true for V/^. To have an approximate comparison to the scaled 
elliptic flow, Fig. IV. 9 shows a scaled hexadecupole flow with Vi/v? vs. Pr/n. The 
is inspired by an observed constant ratio t'4/(t'2)^ (for details see sec. IV. 4.3). 



IV. 2.3.3 Relative strength of both contributions 

While the flow contributions are flxed by the impact parameter, the geometrical contri- 
butions depend on the orientation of the freeze-out ellipse and its eccentricity. To relate 
the initial and the freeze-out eccentricity, I will use a simple linear ansatz 



(IV.14) 



with the parameter Ce < 1. 

To obtain a flrst approximation, let us compare the results from the previous chapter 
with a circular freeze-out (ce = 0) to actual data (see Fig. IV. 27 on page 68). As you 
can see the elliptic flow is generally underpredicted at low pT- Since po only models the 
high-pT damping, it has no effect to change it. The only possibility to raise the pion V2 
high enough would be to lower the quark mass to a few MeV, which not only is very 
questionable, but also fails to reproduce the baryons. 

Hence, to match the data, one seems to need some additional positive contributions 
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Figure IV. 9: Scaled flow coefficients of different liadrons for a circular freeze-out. 
Left: Scaled elliptic flow. Right: Scaled hexadecupole flow. 
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Figure IV. 10: Flow coefficients of light quarks for a circular freeze-out. 
Left: Elliptic flow. Right: Hexadecupole flow. 
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Figure IV. 11: Comparison of the freeze-out eccentricity from eq. (IV. 17) for different Ce 
to STAR calculations [Ada04], where the freeze-out eccentricity has been obtained with 
azimuthally-sensitive HBT and the initial eccentricity was calculated within a Glauber 
model. Uncertainties on the precise nature of space-momentum correlations lead to 30% 
systematic errors on £/. 

from the geometry. This is achieved by setting Cg ~ 0.7. That means the freeze-out 
eccentricity should be about 70% of the initial eccentricity. This also increases the high 
Pt region which is compensated by decreasing po to 0.65. 

This is a first hint that ej is not necessarily zero and needs to be incorporated. This is 
supported by STAR data [Ada04], where the initial and final eccentricity was calculated 
via the Glauber model and HBT interferometry. Their results indicate that the elliptic 
flow cannot quench its initial asymmetry completely. 

As discussed in sec. III. 3.2, the numerical eccentricity e is compared to the glauber 
eccentricity e via 

1 , , 

^ =7:7; 9 respectively (IV. 15) 



Together with the linear scaling of ej = CeCj this relates the initial and the freeze-out 




(IV.16) 
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eccentricity as 



2 



^1 c'ef ^1 2ei + l 



2e, + 1 - 2cJei 



2ei + l 

(IV.17) 



Fig. IV. 11 shows the resuhs of eq. (IV. 17) to the STAR calculations, which gives Ce = 
0.731 ± 0.015 for the best fit case. This good agreement supports the simple linear ansatz 
in eq. (IV. 14) and confirms the first estimate of Cg 0.7. But the large systematic errors 
do not allow a firm conclusion, so £ ^ 0.65 — 0.8 is consistent with errors. 

This data also clearly indicates a positive ej which means that the orientation has not 
changed compared to the initial ellipse and it is still elongated out-of-plane. From what 
I already said in sec. IV. 2.3.1, I conclude that the transverse expansion of the system is 
perpendicular to the surface and drop the surface flow model no. 2 in favor of model 1 
(sec. III. 2.2). 



IV. 2.4 Correlation between r and p 

So far I only studied the elliptic flow with a constant freeze-out time. To be a bit more 
general I assume a linear dependence on the radial coordinate as 

T = Xp + t' and therefore dpT = X. (IV. 18) 

For a given mean freeze-out time tq it follows 

r' = To- ^poX (IV. 19) 

With A = a constant freeze-out time r = tq is recovered. This is a quite simple ansatz, 
but it suffices to study a general dependence. 

The dependence of the yield on A is shown in Fig. IV. 12 for different hadrons. The 
relative deviations rx{dN) = {dNx — dNx=o) /dNx=Q is similiar for mesons and baryons 
and approximately independent of px with about -0.2 for A = 0.3 and 0.2 for A = —0.3. 

The deviations Ax{v2) = ^2^'* — '^2^"'''* of the elliptic flow for different A are shown in 
Fig. IV. 13 for an elliptic freeze-out with c,. = 0.7. For a circular freeze-out the deviations 
are neglectable. So non-zero A has the largest effects on V2 at low px- This behaviour is 
similar for the hexadecupole flow. 
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Figure IV. 12: The relative deviations rx{dN) of the invariant yields for A — ±0.3 to 
A = in central coUisions {b — fm) . 




Figure IV. 13: The deviations A\{v2) of the elliptic flow for A — ±0.3 to A = in mid- 
central collisions {b — 8 fm) for an elliptic freeze-out (ce ~ 0.7). 
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IV. 3 Transverse momentum spectra 

In this section I show hadron production for Au+Au colhsions at \/s = 200 GeV from 
recombination. For the LHC predictions of Pb+Pb colhsions at ^/s = 5.5 TeV I will use 
the same nuclear size po, since the radius only enters an overall factor and there is no data 
yet available to compare to. 



IV. 3.1 Yields 

Fig. IV. 14 shows the transverse momentum spectra of light and strange hadrons compared 
to data from PHENIX [Chu03] and STAR [Abe06, Min08, Abe07, AdaOTa] for different 
centralities. In the mid pt range, the results agree with the data. But below 2 GeV 
and especally below 1 GeV, the data is significantly underpredicted for hadrons with light 
quark content. This is due to the energy violation in the recombination formalism and 
also the contributions from resonance decay are not incorporated. For hadrons with only 
strange quark content, namely (f) and the predictions are in good agreement even down 
to very low px, since for the massive particles the energy violation is not that strong and 
the contributions from resonance decay are minor. The black lines correspond to LHC 
predictions for an impact parameter of 6 = 2 fm. 

Fig. IV. 15 shows the transverse momentum spectrum of charged hadrons compared to 
data from STAR [Ada03] for different values of dp, the parameter for the radial profile of 
the transverse rapidity. While the results are equal at low px-, they differ for joj^ > 3 GeV 
and the best fit is achieved for about dp = 0.5 — 1. But the large deviations from the data 
at low Pt do not allow for a firm conclusion on that parameter. Another analysis is given 
in sec. IV. 4.3.2, where I study the dependence of the flow ratios on this parameter. For 
now I will use dp = 1. The other parameters are Cg = 0.73 and A = —0.3 as discussed 
in the previous sections. The negative A gives a constant factor of about 1.2 while the 
freeze-out eccentricity does not affect the transverse momentum spectra, since it is chosen 
to preserve the total transverse freeze-out area. 



IV. 3.2 Mean transverse momentum 

Fig. IV. 16 shows the mean transverse momentum {px) as a function of the hadron mass 
in the left panel. The results for RHIC (red points) are compared to data from STAR 
[AbeOSa] (blue points). At low hadron mass, namely for pions and kaons, the results for 
the mean pr overestimate the data, since the low pT yield is underestimated as shown in 
the previous paragraph. For protons and heavier hadrons, the deviations for the yield are 
smaller and therefore the predicitions of the mean px are in much better agreement. 

Additionally, in the right panel I show the mean transverse velocity (t'T)L=o = = 

(pt) 

at midrapidity calculated from the mean px- As one can see, this quantity softens 

{mT) 

the deviations for the light mesons. For masses of about 2 GeV the velocity approaches 
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Figure IV. 14: Transverse momentum spectra of hadrons for Au+Au collisions at -y/s — 
200 GeV compared to pion, kaon and proton data [Chu03] from PHENIX (blue points), 
pion and proton data [Abe06], kaon data [Min08], (f) data [Abe07], A and u data [AdaOTa] 
from STAR (red and green points). 
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Figure IV. 15: Comparison of the charged hadron yield for different radial profiles of the 
transverse expansion rapidity (dp) to data from STAR [Ada03]. 



the value of the mean transverse expansion /Jy = 0.55 (black lines) at RHIC. Similiar the 
predictions for the LHC (green points) approach (3t = 0.75. 

The figure also depicts calculations from the equipartition theorem (red/green lines), 
which states 

/ dH \ 

(Prn-Q^) = ^mnkBT (IV.20) 



with the temperature T and the Boltzmann constant ks (which is equal to 1 in natural 
units). Using the energy H = E = p^u^^ from eq. (III. 62) in the local rest frame, this 
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Figure IV. 16: Mean transverse momentum (Left) and velocity (Right) of hadrons (/3) = 
(Pt) /E as a function of the hadron mass at RHIC (red) and LHC (green). Data points 
(blue) are from STAR [AbeOSa] and the black lines are the theoretical expectation within 
the equipartition theorem. 



leads to 



2T = prp — px cos((/? — <P) sinh(?7r) = (pr 



These equations are coupled via E = ^Jj^^^\-j^'-\-rn^ so they need to be solved numeri- 
cally. 

Interestingly, the values for hadron masses above 1 GeV are larger then the results from 
recombination, if one uses the same temperature T = 0.175 GeV and transverse velocity 
Pt = 0.55 or respectively (3t = 0.75. Therefore, the mean velocity (Pt) reaches these 
"limiting" values much later. When using T and Pt as fit parameters, the equipartition 
theorem reproduces the recombination results for T = 0.324 GeV, Pt = 0.0 for RHIC and 
Pt = 0.45 for LHC. 
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Figure IV. 17: Invariant yield ratios for different hadrons as a function of Pt compared to 
data from PHENIX [Adl03] and STAR [Lon04, Min08]. 



IV. 3.3 Hadron ratios 

The large proton/pion ratio at mid px was one of the puzzles that fragmentation failed to 
describe and thus, it is one of the motivations for studying recombination. So, Fig. IV. 17 
shows the ratios of the invariant yields from different hadrons compared to data from 
PHENIX [Adl03] and STAR [Lon04, Min08]. The low pr region is mainly overestimated 
which can be attributed to the failed description of the pT spectra at low pT, due to 
the energy violation. The mid px results are between a mediocre and good agreement, 
apart from the A to ratio which is about a factor of 2 smaller. The main point is 
the prediction of the large proton to pion ratio of about 1. The high px region is badly 
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predicted and would require the combined treatment of recombination and fragmentation 
[Pri03]. 

The black lines correspond to LHC predictions, but with the same strange fugacity 
7s = 0.8 which would have to be extracted from fits and can be expected to be larger. 



IV. 4 Flow coefficients V2 and V4 
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IV. 4 Flow coefficients V2 and 

In this section, I study the eUiptic and hexadecupole flow at RHIC and LHC. I will start 
with the eccentricity dependence of the mean flow coefHcients. The findings will motivate 
the investigation of the flow ratio ^4/(^2)^ which will support the findings in sec. IV. 2.3.3 
about the freeze-out eccentricity. 

This ratio also gives the opportunity to study a r — p correlation and the radial profile of 
the transverse expansion rapidity. The obtained values for the corresponding parameters 
in the sections IV. 4.3.1 and IV. 4.3.2 will be used in the comparison to experimental data 
of the differential flow coefficients. 



IV. 4.1 Mean flow coefficients 

The mean elliptic or hexadecupole flow {vn) at midrapidity is the px integrated average 
folded by the transverse momentum distribution: 

(Vn) = JJ^ (IV.23) 



J dp 



r 



dpi 



IV. 4.2 Eccentricity dependence 

The strength of the flow coefficients Vn depends on the impact parameter or respectively 
the eccentricity. To study the dependence, I will compare the mean, pr-integrated flow 
(vn) to the eccentricity e. Looking at equation eq. (III. 56), one expects the mean V2 to 
scale mainly linear with e and V4 to have a larger quadratic contribution. 

Therefore, the calculations are fitted by i^) = o,e + 6e^ and shown in Figs. IV. 18 
and IV. 19. The fit values support the expected behaviour and can be found in Tab. IV. 1. 



Table IV. 1: Fit values for the function (e) = ae + he^ as shown in Figs. IV. 18 
and IV.19. 
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Figure IV. 18: Mean flow (?;„) as a function of the eccentricity e for mesons (left) and 
baryons (right) from a circular freeze-out (cg = 0). The lines are fltted by (u„) (e) = ae + be^ 
with the values from Tab. IV. 1. 

Upper: Mean elliptic flow. Lower: Mean hexadecupole flow. 
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Figure IV. 19: Mean flow (vn) from an elliptic freeze-out with Cg — 0.73. Tab. IV. 1 shows 
the flt values. 
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IV. 4.2.1 Fluctuations 

For a fixed impact parameter, the eccentricity is expected to fluctuate [\bl08, Bro07]. 
Since the flow components scales with powers of e, the fluctuations in the pT-mtegrated 
{vn) can be related to the eccentricity fluctuations. To study these fluctuations, I take the 
fits from table IV. 1, where I used 

{vn){e)=ae + he'^. (IV.24) 

I assume a gaussian distribution of the eccentricity around the mean expected eccentricity 
£ as 



pie) = iVexp 



2 a2(e) 



(IV.25) 



with a width a that generally can depend on the mean eccentricity. The mean of the 
elliptic flow distribution can then be calculated as 

1^ = b{a'^{e)+i^)+ae (IV.26) 
and the fluctuations are 

a{{vn)) = a(£)^262(2(^-6a2(£))+a2 (IV.27) 

where (x) is the pT-integrated mean and x denotes the fluctuation mean. 
The relative fluctuations for a purely linear dependence {b = 0) simplify to 



(IV.28) 



For a purely quadratic dependence (a = 0) the relative fluctuations scale with Vb as 

^^^""^^~2v/5f#-^(^) 1 (IV.29) 




{Vn) V 

Looking at the fit values in table IV. 1, the quadratic contributions for the elliptic flow 
are small compared to the linear ones (a > b) and vice versa for the hexadecupole flow 
(6> a). 

The relative fluctuations of the eccentricity are predicted to be approximately constant 

with ^ ^ 0.4 [Dre07, Alv07, Sor07]. As can be seen in Figs. IV.20 and IV.21, the 
relative fluctuations of the elliptic flow is almost independent of the particle type and also 
comparable to the fluctuations of the eccentricity as eq. IV.28 suggests. On the other 
hand the relative hexadecupole fluctuations are much smaller and approch the value from 
eq. IV.29 for large e. The divergence of cr{v4)/{v4) at low px is due to the negative fit 



60 



IV Results and predictions from recombination 




value of a, so that the denominator becomes zero at a finite e. 
IV. 4.3 Flow ratio 

As was shown in the last section, for a circular freeze-out the different flow coefficients Vn 
will scale with the eccentricity mainly as e"/^. Therefore ratios of these coefficients are 
an interesting probe to study, since the flow contributions will be insensitive to the initial 
geometry. The geometrical contributions from an elliptic freeze-out on the other hand can 
have great influence on these ratios. So they offer the possibility to verify the statements 
about the freeze-out eccentricity. One can also take a closer look on the influence of a 
T — p correlation (A / 0) and the radial profile of transverse expansion (dp). 

I will focus on the ratio — ^, that I already discussed shortly in framework of the CQNS 
(sec. III. 4.3.2). For a quick approximation of the value for the quarks I compare the 
prefactors from the transverse rapidity (eq. (III. 56)): f 2 ~ e and ~ -e^, which yields 
-1 quark ^ 

-. This will turn out to be too small. Inserting a more realistic value of 1 
vj\ 4 

in the scaling law from section III. 4.3.2 gives constant values of 



V4 



meson ^ 



4 + 2 = 4 (^^-^^^ 



IV. 4 Flow coefficients V2 and V4 



61 





0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 
Pj [GeV] P-p [GeV] 

Figure IV. 22: The ratio ^4/(^2)^ for charged hadrons from an elhpsoidal freeze-out with 
different Ce and two impact parameters. 
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and 



-\ barvon 



1 1 _ 2 
3 ^ 3 ~ 3 



(IV.31) 




So the simple cqns formulas suggest a ratio around 0.7. But these values are much smaller 
than the experimental values from STAR [Bai07] (with ~ 1.2) and PHENIX [Hua08] (with 
~ 0.9 ). Also the full calculations give a similiar value of about 0.7 from mid to high px- 
That indicates that we have to consider contributions from geometrical effects. 

Fig. IV. 22 shows the flow ratio for charged hadrons. The lines correspond to different 
percentages Cg of the freeze-out eccentricity to the initial one. The left figure corresponds 
to an impact parameter of 6 = 2 fm, the right one to 6 = 8 fm. The good agreement 
between both supports the predicted independence from the eccentricity. 

At low pt the Vi/{v2)^ ratio is much greater than 1 for all percentages of the freeze-out 
eccentricity. At mid pT the effect of an elliptic freeze-out is most clearly visible. While 
for Cg = the ratio is about 0.6, it rises to about 1 for Cg = 0.73. The three curves with 
Cg = 0.65,0.73 and 0.8, which are the estimates from sec. IV. 2.3.3 within errors, only differ 
at low Pt- For these three values the ratio is in the range of the experimental data. So 
to compare the results to experimental data, I will take the best fit case of Cg = 0.73 and 
in the next two paragraphs I look at the dependence on two other parameters, namely A 
and dp. 

IV. 4.3.1 Influence of a r — p correlation 

This is dependence is studied by varying the parameter A (see section IV. 2.4). Therefore 
I fix dp = 1, which means a linear growth if the transverse rapidity. The upper panel in 
Fig. IV. 23 shows the ratio for charged hadrons at 6 = 2 fm for A = 0, ± 0.3 compared 
to data from STAR [Bai07, Pos04]. A constant freeze-out time, i.e A = 0, seems to give 
the best agreement. This different to the lower panel with t'i/(i'2)^ at = 8 fm for pions, 
protons and kaons, where at high pT the value is independent of A, but the predictions for 
A = at low Pt overestimate the data. A very good agreement at least for the mesons can 
be found with A = —0.3. A (negatively) larger A for baryons could increase the agreement 
at low Pt, but the very sharp drop of the proton ratio can not be accounted for. And it 
also would enhance the discrepancy to the charged hadron data from STAR. 

A negative A describes a fireball, where the outer particles freeze-out earlier (tq than the 
inner ones (rj. With A = —0.3 and a mean freeze-out time tq = 5 fm, this corresponds 
to To = 4.1 fm and Tj = 6.8 fm. The data does not allow to draw a firm conclusion about 
the value of A, since it also could (or even should) depend on the particle species or on 
their in-medium cross sections respectively. But the case of a large positive A, where the 
To > Ti, can be discarded. 



IV. 4 Flow coefficients V2 and V4 



63 



> 




J TJ I I I L 



> 



1.4 
1.2 

1 

0.8 
0.6 



\ 1 M 1^1 1 1 




\ \ 




\ \ 












1 : 





1.5 2 2.5 3 3.5 
Pj[GeV] 



X = 0.0 

A, = -0.3 b = 8fm 
A, = 0.3 

PHENIX data, cent. 20-60% 



0.5 1 1.5 2 2.5 3 3.5 
Pt [GeV] 

Figure IV. 23: The flow ratio W4/(w2)^ for an ellipsoidal freeze-out with Ce = 0.73. The 
influence of a r — p correlation is depicted with A = 0, ± 0.2. Upper: Ratio at & = 2 fm for 
charged hadrons compared to data from STAR [Bai07, Pos04]. Lower: Ratio at 5 = 8 fm for 
pions, kaons and protons compared to data from PHENIX [Hua08]. 
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IV. 4.3.2 Influence of the radial rapidity profile 

The radial profile of the transverse rapidity is set by the parameter dp (see section IV. 2.2.4). 
So this time I use a constant A = 0. 

The behaviour for different dp is similiar to the previous paragraph. At high pT all 
three curves give the same result and only at low pr the values are modified. The data on 
charged hadrons (Fig. IV. 24, upper panel) is best reproduced with dp = 1 (and A = 0) as 
in the previous paragraph. The data from PHENIX on pions and kaons seems to favour 
a smaller dp of about 0.5, while proton data at low pT points to dp > 2, where dp = 2 still 
gives a large overestimation. 

So generally both parameters can be used to fine tune the V4^/{v2f' ratio and this would 
preclude a determination of these parameters from this observable. But there are two 
reason to fix dp at a value of 1 and use A as a free parameter: 

• The low pt pion data could be explained by a lower value of A for baryons. A 
larger value of dp on the other hand is not only unreasonable, but also largely fails 
to predict the meson data. And since dp is a blast-wave parameter is would be 
inconsistent to let it depend on the particle species. 

• The outer particles in the shell of the fireball suffer less collisions, since they can 
expand into the vacuum, and therefore one expects them to freeze-out earlier than 
the particles in the inner fireball. From that, one would expect a negative A < 0, 
which is supported by hydrodynamical calculations [Kis06, Cho06]. 

IV. 4.4 DifFerential flow 

In this section I will show differential elliptic and hexadecupole flow Vnipr) compared to 
data from RHIC. I will use the parameters as extracted from the previous sections as 
Ce = 0.73, dp = I and A = —0.3 if not specified otherwise. To show the relevance of 
these parameters. Fig. IV. 25 compares the elliptic flow of charged hadrons with an impact 
parameter of 6 = 8 fm for three different freeze-out scenarios: circular freeze-out with 
Ce = 0, A = (blue line), elliptic freeze-out with Cg = 0.73, A = (green line) and elliptic 
freeze-out with c,. = 0.73, A = —0.3 (red line). Again, a circular freeze-out can be safely 
discarded. An elliptic freeze-out with a radial decrease of the freeze-out time (A < 0) is the 
best description of the STAR data [Tan04] and only gives a slight underestimation at low 
Pt- The data points correspond to the event plane (squares) and two-particle cumulant 
method (circles) for the determination of the elliptic flow. 

Fig. IV. 26 compares the differential elliptic flow V2{pt) of identified hadrons to data from 
RHIC from different centralities. The red data points are from STAR with a centrality of 
40-80% and the green ones from PHENIX with 20-60% centrality. They are compared to 
an impact parameter of 9 fm (red lines) and 8 fm (green lines) respectively. 

The agreement is quite good. While the behaviour at high px changes with the phe- 
nomenological parameter pQ, the precise predictions at mid and low pr confirm the rele- 
vance of the parameter Cg. To show that a large freeze-out eccentricity of about 70% of 
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Figure IV. 24: The flow ratio W4/(w2)^ for an ellipsoidal freeze-out with Ce — 0.73. The 
behaviour for different radial profiles of the transverse rapidity is depicted with dp = 0.5,1 
and 2. Upper: Ratio at 6 = 2 fm for charged hadrons compared to data from STAR [Bai07, 
Pos04]. Lower: Ratio at 6 = 8 fm for pions, kaons and protons compared to data from 
PHENIX [Hua08]. 
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Figure IV. 25: Elliptic flow of charged hadrons as a function of pt for Ce — (blue line) 
and Ce — 0.73 (red line (A = —0.3), green line (A = 0)) for an impact parameter 6 = 8 fm 
compared to data from STAR, cent. 20-60% [Tan04] (points). 

the initial one is not only need for the flow ratio from the previous section, but also for 
the differential elliptic flow, Fig. IV. 27 gives the same comparison for Ce = and A = 0. 

So for a circular freeze-out the low and mid px range is generally underpredicted for 
all hadrons and the high pT region only agrees, because of a larger po. As can be seen, a 
remaining freeze-out eccentricity can not be neglected. 

The predictions for LHC (black lines) are also shown. A first striking observation is 
that the values are generally smaller. A similar pattern was also observed within a parton 
transport approach, if the viscosity was set to the ADS/CFT limit [Mol07]. It seems as 
the elliptic flow reaches a maximum somewhere between these two center of mass energies 
which is due to the increased transverse expansion velocity from 0.55 (RHIC) to (3t = 0.75 
(LHC). I have already shown this strong dependence on Pt in section IV. 2.2.3. At low px 
the elliptic flow becomes even negative for (multi-) strange particles which will be exploited 
in the following sections IV. 4.6 and IV. 4.8. A discussion about the negative flow then 
follows in section IV. 4.9. 

For the fourth fourier coefficient v^, experimental data is rare up to now. Fig. IV. 28 
shows the hexadecupole flow of different hadrons for an elliptic freeze-out compared to 
the available data from PHENIX [Hua08]. The pion results agree very well with the data 
while the kaon and proton data is slightly overestimated at low pT- Also the predictions 
for the LHC (blue lines) are shown. 
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Figure IV. 26: Elliptic flow of identifled hadrons as a function of pt with Cg — 0.73 and 
A = —0.3 for impact parameters of 9 fm (red lines) and 8 fm (green lines) compared to 
data from STAR, cent. 40-80% [Abe08b] (red points) and PHENIX, cent. 20-60% [Afa07] 
(green points). The black Hnes correspond to the predictions for LHC {^/s = 5.5 TeV) with 
6 = 8 fm. 



IV. 4.5 Deuteron flow 

Due to their small binding energy, deuterons are expected not to survive the early phase, 
but to be produced by proton-neutron coalescence at freeze-out [Sch97, Nys04]. Thus, the 
invariant yield can be written as 



Ed-^ = B2 ( Ep-^ ) (IV.32) 
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Figure IV. 27: Elliptic flow as a function of pt /or a circular freeze- out with Ce — 0.0 
and A = for impact parameters of 9 fm (red lines) and 8 fm (green lines) compared to 
data from STAR, cent. 40-80% [Abe08b] (red points) and PHENIX, cent. 20-60% [Afa07] 
(green points). The black Hnes correspond to the predictions for LHC (i/s = 5.5 TeV) with 
6 = 8 fm. 



with pd = 2pp and a coalescence parameter B2 which will depend on the freeze-out volume. 
Hence, the elliptic flow is related by 

.f(P^)=2.?)(^). (IV.33) 

But by inserting the CQNS for the proton elliptic flow, the expected scaling of the deuteron 
flow with the quark flow is obtained: 

.f(P^)=64^)(f). (IV.34) 
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Figure IV. 28: Hexadecupole flow as a function of pt with Ce = 0.73 and A — —0.3 for 
an impact parameter of 8 fm (lines) compared to data from PHENIX, cent. 20-60% [Hua08] 
(points). The blue lines correspond to the predictions for LHC (-/s = 5.5 TeV) with b = 
8 fm. 



Thus, the general behaviour of the deuteron elhptic flow seems to be independent of the 
production process: uud+udd — > p+n d (after the freeze-out of p and n) or uuuddd d 
(at hadronization). 

Fig. IV. 29 compares both cases with realistic calculations at 6 = 8 fm to data from 
PHENIX [Afa07] (cent. 20-60%) and STAR [Liu07] (minimum bias). The red line corre- 
sponds to the scenario of 6 recombining quarks at hadronization which agrees very well 
with the data. Only at very low px, it does not follow the STAR data of a slightly nega- 
tive V2- The green line depicts the scaled proton elliptic flow 2v^^ {Pt/2) which does not 
give the same results but follows the same trend as expected. It overestimates the data 
which may be due to the "delta-shaped wavefunction" approximation Pp = Pn = Pd/2. So, 
while the study of the transverse momentum spectra requires a dynamical treatment with 
hadronic degrees of freedom, the elliptic flow can be described by quark recombination 
and thus follows the CQNS. Additionally, the black line shows the prediction for the LHC 
for the quark recombination scenario. 
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Figure IV. 29: Elliptic flow of deuterons from 6 recombining quarks (red line) or scaled 
proton V2 (green line) compared to PHENIX [Afa07] and STAR [Liu07] data. The black 
lines correspond to the LHC prediction. 



IV. 4.6 Heavy quark flow 

A large elliptic flow implies, from a hydrodynamical point of view, a rapid thermalization 
of the fireball and a strong collective flow created at the QGP stage. This assumption 
is well established and the large flow of light quarks is reflected in the good agreement 
with the V2 data. Heavy quarks on the other side have a much larger mass then the light 
quarks and therefore it is an open question, if charm or even bottom quarks participate 
in the collective expansion. 

The recently published data on J/ip from PHENIX offer the possibility to investigate 
whether also the charm quark does locally equilibrate and therefore follows the flow of the 
light quarks. The most prominent feature of this preliminary data, despite its large errors, 
is the negative V2 value at px = 1-5 GeV. Within the CQNS this implies also a negative 
elliptic flow for the charm quark. Since the multistrange f2 also shows a slightly negative 
V2 for a mean transverse expansion velocity of /? = 0.55c, the much heavier charm quark 
can be expected to show an even more pronounced behaviour, when assuming a similar 
transverse expansion. 

The elliptic flow for J/ip is depicted in Fig. IV. 30 together with the data from PHENIX 
[Sil08]. The left panel is calculated with the parameter A = and shows different expansion 
velocities and the right panel compares /? = 0.55c for different values of A. 

Looking at the left panel, the lines from top to bottom correspond to (3 = 0.3, 0.4, 0.5 
and 0.6c with Ce = 0.73, while the last line (black) is for (3 = 0.55 from a circular freeze- 
out (ce = 0). So a large charm flow that is equal to the light quark flow seems to 
produce large out-of plane elliptic flow {v2 < 0, black line), but this is compensated by 
the positive geometrical contributions from the elliptic freeze-out. Nevertheless, small 
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Figure IV. 30: Estimating the charm flow from the eUiptic flow of J/tp compared to data 
from PHENIX [Sil08]. 

Left: Comparison for increasing mean expansion velocities f3 from top to bottom with Ce = 
0.73 and A = 0. The lowest Hne (black) corresponds to /5 = 0.55c and Ce = 0. 
Right: Comparison for different A with /? = 0.55c and Ce — 0.73. 



transverse expansion velocities seem to be discarded due to the large positive V2 at pT = 
1.5 GeV, while a thermalized charm with /? « 0.5 — 0.6c is at least consistent with zero 
flow at pt = 1-5 GeV. I would like to mention that the discussion from [Kri07] is for a 
circular freeze-out only, since we had employed only the simpified equations from [Fri03]. 

Prom the flow ratio section I predicted a negative A. This delayed freeze-out of the 
particles in the inner fireball can be explained by in-medium rescattering. But the J/ip 
is expected to have a low cross-section, i.e. it can escape the system quite undisturbed. 
So the J/ip will freeze-out early, while the system is still expanding. Therefore one could 
assume a zero or even positive A. This is depicted in the right panel of Fig. IV. 30 for a 
mean expansion velocity /3 = 0.55c, where the data clearly seems to favor a positive A 
(blue line). 

The large errorbars do not allow to draw any firm conclusions, but I want to state that 
the data is consistent within errors with a charm quark flow of /3 = 0.55c equal to the light 
quarks and a positive A « 0.6 (r^ < Tq). At least one can extract a lower bound on the 
expansion velocity of about 0.4c. So if more precise data will still support the negative V2, 
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Figure IV. 31: Electron elliptic flow data from heavy flavor decay from PHENIX [AdaOTb] 
compared to Dq V2 for different charm transverse velocities (3c = 0.3c, 0.4c, 0.5c (Hnes, top to 
bottom) with Ce = 0.73 and A = 0. 



I conclude from this observation that charm quarks reach a substantial amount of local 
kinetic equilibration. 

Another possibility is to study D-meson elliptic flow. There, the positive light quark 
V2 will compete with the negative one from the charm quark. So far, there is no data on 
D-meson elliptic flow available, but in the near future the Heavy Flavor Tracker (HFT) 
from STAR will close this gap. Preliminary, one could compare to V2 data of electrons 
from heavy flavor decay. This is no direct probe of the charm flow, since there will be 
contributions from -B-mesons and the decay kinematics might smear out the charm flow 
signal. However, calculations within transport theory [Gre04] predict the Z)-meson elliptic 
flow to be similiar to the heavy-electron V2. 

Fig. IV. 31 compares the electron elliptic flow data from PHENIX [AdaOTb] to Dq for 
different charm transverse velocities (3c = 0.3c,0.4c,0.5c (lines, top to bottom) with Ce = 
0.73 and A = 0. The Dq elliptic flow agrees with the electron V2 data at low pT, but 
does not allow for a distinction between the different transverse velocities. While the data 
follows the curve for a /3t ~ 0.3 below px = I GeV, above 1 GeV it is better predicted by 
a large expansion velocity of (3t ~ 0.5. But at high pT, the V2 of the electrons is generally 
about a factor of 2 smaller then that of the Dq. That might be due to an early onset of 
the fragmentation regime for the charm. Again, there is no clear conclusion, but the data 
does not seem to contradict a large charm transverse velocity, which is comparable to the 
light quark velocity of Pt = 0.55c, as far as the heavy electron V2 can be compared to the 
Dq elliptic flow. 

As a finally comparison. Fig. IV. 32 shows the elliptic flow for charm and bottom mesons 
at RHIC (left) and LHC (right) for Ce = 0.73 and A = 0.3. 
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Figure IV. 32: Elliptic flow of heavy mesons with charm and bottom quark content at 
RHIC (left) and LHC (right) for & = 8 fm, = 0.73 and A = 0.3. 

IV. 4.7 Centrality dependence 

Fig. IV. 33 shows the mean elliptic flow of charged hadrons as a function of the centrality. 
The left panel is for an elliptic freeze-out and the right for a circular one. Although the 
differential V2 shows very good agreement with the data as shown above, the results for 
(^2) (left panel, red line) generally overestimate the data from STAR [AbeOSb]. This is 
due to the yield at low pT, which is to small because of energy violation and the absence of 
resonance feed-down as discussed in section IV. 3.1. Therefore the mean V2 weighted with 
the yield favors larger V2 values. I account for that by scaling the results with a constant 
factor of 0.75 (green line), so that they are in agreement with the data for the central and 
mid-central collisions, but for peripheral collisions the data is still overestimated. This is 
due to the used eccentricity 6/(4i2^) (eq. (III. 54)), which follows the glauber eccentricity 
for central and mid-central collisions, but gives wrong results for peripheral collisions 
(Fig. III.2). 

The mean V2 from a circular freeze-out fails to describe the data and underestimates it 
as expected. Only to have a comparison between both cases, I will scale the mean elliptic 
flow value for Cg = by a constant factor of 1.3 (right panel, green line) so that it fits 
the data for central to mid-central collisions. In the next section I will use both factors 
implicitly. 

Due to the lack of similiar data for the hexadecupole flow, I compare the results to 
V4 of charged hadrons for different centrality classes from PHENIX [Iss08] in Fig. IV. 34. 
The agreement with the data is quite good within the uncertainty of relation between 
the impact parameter and the centrality. But an underestimation of the mean (^4) for 
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Figure IV. 33: Mean elliptic flow for charged hadrons as a function of the centrality com- 
pared to data from STAR [AbeOSb]. 

Left: {V2) from an elliptic freeze-out. Right: {V2) from a circular freeze-out. 



peripheral collisions can be expect similiar to the mean elliptic flow. 
IV. 4.8 dependence 

The effect of a negative elliptic flow at low px at LHC will be visible even for hadrons with 
light quark content, like protons as shown in Fig. IV. 26. Therefore, it is interesting to 
look at the mean elliptic flow as a function of the center of mass energy ^/s. To compare 
the results from an elliptic and circular freeze-out. Fig. IV. 36 uses the correction factors 
from the previous section. 

The shown value for LHC in Fig. IV. 36 is just a linear interpolation, but the general 
expectation in the heavy ion community is an monotonic increase from RHIC energies to 
LHC and beyond, or at least a saturation at some finite value. In this framework, the 
value of (^2) depends on the interplay of the strength of the negative V2 and the increasing 
mean pT, since (^2) ~ V2 ((pt))- But as can be seen in Fig. IV. 36, the increase in the 
mean pT does not seem to compensate the increasing effect of the negative V2 at low pT- 
Depending on the used freeze-out eccentricity, the maximum of the mean elliptic flow is 
reached somewhere between RHIC and LHC energies, while a larger eccentricity (a larger 
Ce) delays the extremal point to higher cm. energies. For the best fit case of Cg = 0.73 
the value at LHC is maximal and then starts to drop; for a circular freeze-out with Cg = 
it would have already reached its maximum at RHIC. This may mean that even if this 
striking prediction is correct, there will be no clear sign of a decreasing mean elliptic flow 
for charged hadrons at LHC. But the mean elliptic and hexadecupole flow of identified 
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Figure IV. 34: Hexadecupole flow as a function of with Ce = 0.73 and A = —0.3 for im- 
pact parameter of 6,7 and 8 fm (lines) compared to data from PHENIX (points) of different 
centrality classes [Iss08]. 

hadrons can be a much better probe. As can be seen in Fig. IV. 35, the V2 and U4 of pions, 
(/>'s and protons reach their maximum at LHC energies. And while the V2 of the J/^p and 
f2 only indicate a small decrease, the V4 of these heavy hadrons will show a visible drop 
when going from ^ = 200 GeV (RHIC) to 5.5 TeV (LHC). 

A rather critical assumption in this context is the applicability of the recombination 
approach for the elliptic flow for small transverse momenta on the order of px < 1 GeV. I 
want to emphasise that the result of the decreasing mean elliptic flow (^2) at LHC is not 
affected by the validity of this assumption, because {V2) « V2 {{pt)) and (pt) > 1 GeV 
in LHC regime. However, to show the robustness of the prediction Fig. IV. 37 depicts 
the elliptic flow at a fixed px as a function of ^/s with an impact parameter 6 = 6 fm. 
With Pt = 0.6, 1 and 2 GeV the V2 exhibits the same drop as the mean elliptic flow from 
Fig. IV. 36. The data points are taken from [Adl05] with a centrality of 13-26%. They 
confirm the observation of an elliptic flow saturation at ^/s = 200 GeV. 

As can be seen in Fig. IV. 38, the mean hexadecupole flow shows a similar dependence 
on the center of mass energy. In the case of an elliptic freeze-out (red line), the increase 
from RHIC to LHC is much more pronounced as for the elliptic flow, but also the drop 
beyond its maximum at LHC is much stronger. In the case of a circular freeze-out (green 
line), the predicted value is already negative at LHC. 

Since both flow coefficients show a similiar behaviour as a function of I will also take 

a look at the ratio -r-^L. Fig. IV. 39 shows this ratio for charged hadrons as a function of 

{V2)^ 

the center of mass energy at an impact parameter of 6 = 1 fm (solid lines) and 6 = 6 fm 
(dashed lines). The results for both impact parameters are nearly equal, hence this mean 
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Figure IV. 35: Comparison of (1^2) and (■i;4) for different hadrons from -^s = 100 GeV to 
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Figure IV. 36: Comparison of (^2) for charged hadrons as a function of center of mass 
energy for Au+Au/Pb+Pb reactions at an impact parameter of 6 = 6 fm. The calcula- 
tions are scaled with the constant factors from Fig. IV. 33 as discussed in sec. IV. 4.7. Data 
points and the extrapolation to LHC are taken from [Bor08]. 
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Figure IV. 37: Elliptic flow 112 of charged hadrons at 6 = 6 and flxed pr as a function of 
center of mass energy for pr = 0.6, 1 and 2 GeV compared to data from PHENIX [Adl05]. 
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Figure IV. 38: Comparison of (U4) for charged hadrons as a function of center of mass 
energy for Au+Au/Pb+Pb reactions at an impact parameter of & = 6 fm. 
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ratio is also independent of the centrality, like the pT dependent ratio in section IV. 4.3. 
For a circular freeze-out (green lines), the ratio is approximately constant between 1.3-1.5 
up to RHIC energies and is steeply dropping when going to LHC while for an elliptic 
freeze-out (red lines), it is constant at about 1.3 up to LHC energies. 

IV. 4.9 Analyzing the negative elliptic flow 

At low center of mass energies, the elliptic flow is generally negative, since the spectators 
in a non-central collision are blocking the in-plane flow. This behaviour, also known 
as anti-flow, is not observed at high bombarding energies, since the nuclei are strongly 
lorentz-contr acted. Then the flow asymmetry only depend on the spatial asymmetry as 
discussed. 

The predicted negative elliptic flow at low transverse momenta is different to anti-flow. 
It has also been found in previous exploratory studies and seem to be a general feature of 
the blast-wave like flow profile at high transverse velocities [Vol97, HuoOl, Vol02, Ret04, 
Pra05]. One might argue that this is an artefact of the blast-wave peak and will not 
survive in more realistic calculations, however, also transport simulations indicate slightly 
negative V2 values for heavy particles at very low pT [Ble02] and the AMPT model predicts 
a negative V2 at least for heavy charm and bottom quarks [Ko07]. Thus, the qualitative 
behaviour of a negative V2 is a well-known observation. The magnitude, however, of this 
effect and the particle species affected by it depend on the mass, the amount of transverse 
flow and the decoupling hyper-surface of this individual particle species. So the surprising 
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observation of a drop in the mean elliptic flow at LHC and an even negative value beyond 
is not the existence of the elliptic 'anti-flow' but the quantitative strength and influence 
on the light quark sector of this effect at LHC. 

To obtain an analytic expression which explains the negative V2 is to consider only the 
in-plane {(p = 0) and out-of-plane = 7r/2) directions (similar to the analysis performed 
in [HuoOl]). The (f) integration breaks down (eq. (III. 80)) and the elliptic flow is then 
given by 

For pt ^ the argument of the Bessel functions In (see eq. (III. 70)) goes to zero and In 
becomes constant with J2 — > and Iq ^ 1. Therefore they are independent of the angle. 
This leads to 

hm vH ) ^2 [mcosh(r?° (1 + e))/T] - [mcosh(7?° (1 - e))/T] 
p ™o ^2 ^P^> 7o Ki [m cosh(r/0 (1 + e))/T] + Ki [m cos\x{r{^{l - e))/T] ^ ' ' 

Since Ki is a monotonically decreasing function, the numerator and thus the elliptic flow is 
negative (for some small transverse momenta). The specific values depend on the mass, the 
mean flow rapidity r/^^ and the temperature. For increasing mass (e.g. charmed mesons) 
or increasing transverse flow rapidity (e.g. LHC and beyond) or decreasing temperature, 
the elliptic flow will become more negative. 

Another way to understand this effect in more detail one could look directly at the 
thermal quark-spectrum with the energy (eq. (III. 62)) 

E = rriT cosh(y — 77) cosh rjT — pr 003(99 — (j)) sinh r]T (IV. 37) 

For simplicity let us look at midrapidity (y = 0) and for high q the spectrum is very low, 
so I consider the region around 7] = 0. For high c.m.-energies, when the source is highly 
boosted transversally, the particles will mainly be emitted in the direction in which the 
fireball flies, so one can simplify even more and set (p = 4>. Because in this case there is 
no longitudinal momentum I can replace the momentum with rj^, the transverse rapidity 
of the parton, as 

Pt =m sinh r)j^ and (IV. 38) 

rriT =\J Tnn? + p'^ = m\Jl + sinh^ ry^, 

=mcosh77f. (IV.39) 

So the energy of the quark in the transverse moving source is given by 



E = mcosh [rij'{l -|- £COs(2(^)) — ryf,] 



(IV.40) 
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IV Results and predictions from recombination 



For a fixed quark rapidity rjj, < rfj (low pt) the energy of the quarks emitted in-plane is 
higher than the energy of quarks emitted out-of-plane. With the thermal spectrum more 
energy means less particles, therefore a negative V2. At rij. = ■qj, one would expect the 
zero-crossing, and above a positive V2. 



V Conclusion 



When studying the elementary building blocks of our universe throughout history, higher 
and higher energies were needed to resolve the ever smaller structures. Now, physicists 
are searching for the so-called Quark Gluon Plasma (QGP), a state of matter where the 
quarks are not confined inside the hadrons but free particles. The transition to this new 
phase of the strongly interaction matter is probed within Heavy Ion Collisions (HIC) at the 
Relativistic Heavy Ion Collider (RHIC) and soon with the Large Hadron Collider (LHC). 

One major problem in this context is the question if we have created a QGP in a collision, 
since the hadronization happens on very small timescales, so that only the hadrons from 
freeze-out can be measured. But also a dynamical description within QCD, of the collision 
in general and the hadronization process in particular, is not expected to be achieved in 
the near future. Thus, we depend on phenomenological models. 

The description via pQCD with the use of measured, parameterized fragmentation func- 
tions is very successful for proton-proton collisions, but in HIC with a dense phase space 
this approach is only applicable at high pT- With the first results from RHIC which pQCD 
failed to describe, a phenomenological model named recombination became popular. 

One of the most promising observables to study the creation of a QGP is the elliptic flow 
which is the second fourier coefficient V2 of the invariant yield. It measures the azimuthal 
asymmetry of the transverse momentum which orginates from the spatial asymmetry in 
non-central collisions. Since it shows a self- quenching behaviour, the measured elliptic 
flow of the hadrons is mainly sensitive to the initial stage of the collision. Therefore, it 
would depend on the partonic elliptic flow from the QGP. 

The success of recombination is mainly based on the prediction of a universal scaling 
law of the elliptic flow which connects the hadron elliptic flow directly with the quark 
elliptic flow. This predicted constituent quark number scaling (CQNS) is experimentally 
well observed and the results can describe the data very accurately at low and mid pr- 
While the qualitative behaviour of the particle species dependence of V2 follows the CQNS 
law of recombination, the quantitative results depend on the two additional "ingredients": 
the quark density distribution and the hadronization/freeze-out hypersurface. 

In this thesis I have shown that these two additionial inputs have strong influence on the 
flow coefflcients. For the quark density distribution I employed the blast-wave model which 
is inspired by hydrodynamics. It describes a transversally expanding, locally thermalised 
fireball. A simple geometrical argument to parameterize the pressure gradient suffices to 
describe not only the elliptic fiow t'2 but also the next coefficient Vi. A result of the blast- 
wave flow proflle is the observed mass scaling of the elliptic flow and especally a negative 
flow coefficient at low pT for massive particles and/or high center of mass energies. Such a 
negative value can also been seen in the preliminary J/ip elliptic flow data from PHENIX 
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V Conclusion 



which seems to support the predictions. The most striking prediction of recombination in 
this framework using the blast-wave model is then a dropping mean V2 at or beyond LHC 
energies. 

The hadronization occurs within a 3-dimensional, time-dependent volume and coincides 
with the kinematic freeze-out, since there is no hadronic scattering phase. The detailed 
study of this freeze-out hypersurface would require a dynamical treatment. But to identify 
a general influence on the spectra, this analytical approach is sufficient, since also in 
hydrodynamics the process of particle freeze-out is a highly non-trivial problem. As I 
have shown, a remaining spatial eccentricity at freeze-out and a radial dependence of the 
freeze-out time have the greatest impact on the flow coefficients, and especally on their 
ratio 1^4/(1^2)^- While a simple circular freeze-out underestimates this flow ratio by a factor 
of 2, an elliptic freeze-out introduces additional contributions to both flow coefficients and 
can describe the experimental data very well. The behaviour of the flow ratio at low pT 
is then affected by the radial dependence of the freeze-out time. 

The presented results in this thesis indicate that recombination is the dominant mech- 
anism for hadron production at mid pr while the flow coefficients and their constituent 
quark number scaling are a strong sign for an early partonic stage of the collision. Despite 
its phenomenological nature, this analytical approach to recombination as the hadroniza- 
tion mechanism in HIC, together with an appropriate quark density and freeze-out hy- 
persurface, describes not only the qualitative constituent quark number scaling of the 
experimental data, but also the quantitative mass scaling of the elliptic and hexadecupole 
flow and the large ratio of these coefficients with great detail. Based on this good agree- 
ments, the striking predicitions for the LHC should be considered in more detail within a 
dynamical model. 
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A Analytical derivations 



A. 1 MIT bag model 

To derive an analytic expression for the phase boundary within the MIT bag model, I 
assume an ideal, relativistic gas of quarks and gluons: 

mq = nig = \p\ = p = E (A.l) 

Thus, the pressure P = -e is related to the energy density e, which can be calculated as 



poo 

£=g E- n{E) dE 
Jo 



I 



r.. \ — n-r (A.2) 



(27r) 

47r f°° 3 1 

" (2^ Jo ^ exp [P{p - BifiB] ± 1 

with the degeneracy factor g, the inverse temperature /3 = 1/T, the baryonnumber Bi. 
The (+) sign is for quarks (fermions) and the (— ) sign for gluons (bosons). Using the 
substitution 

x=p{p-pBfiB) ^P= - + BfiB = T{x + pBiiB), -£ = T (A.3) 

one obtains 

For the gluons g with Bg = 0, it follows 
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A Analytical derivations 



For the quarks q und antiquarks q with Bg = —Bg = 1/3 follows 



~l" — 9q 



--9q 



27r2 
2^ 
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Jo + 1 

^ {X+ PBgfiB)^ 

J-l: 
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J+l 



{X - f3Bq^iBf 



^ {x + (3BqHBf + {x-pBgllBf 



+ 1 



dx 



ax— / ^ da; 



+ 1 



(A.6) 



In the first term, the brakets are expanded and in the third term, is pulled out of the 
denominator and the substitution x — ^ —x is applied: 



Eg + £q —Qq 



°° 2x3 + GxipBgHB? 
+ 1 
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+ 1 
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e^dx 
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120 ' 4 ' 87r2 
The phase transition occurs when the pressure P is equal to the bag pressure B: 

B=P = l{eq + eg + eg) 



(A.7) 



0=T^.C + T'. {2{BgHBfgq) + 3-^^^^ 



- 245 
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-iBgllB)^gg + J iBgllB)% (^ff^ " ^) + ^^^^ 



(A.8) 



TT 



with C = — {7gq + 4(7g). The bag pressure B is then fixed by defining the critical tem- 
15 

perature Tq at zero chemical potential as 

Tc:=T = (/^ (A.9) 



MB=0 



C 



A. 2 Integrals of the flow coefficients 



91 



which I choose to be Tq = 175 MeV. Therefore, the critical temperature as a function of 
the baryo-chemical potential can be written as 



B 



\ 



1 

c 



(A.IO) 



A. 2 Integrals of the flow coefficients 

Here, I will show how to solve the integrals of the fourier coefficients from eq. III. 65: 
Vn =^ j d(?!)COs(r 



dAT 



--C I dp I dip I dr]T{p)p 

Jo Jo J-oo 



dT 



f\ip)mT cosh(y - r?) - —pT (/ cos(/3' - 4>) + f sin(/3' - ct>)) 
exp [— {niT cosh(y — ry) cosh -qx — Pt cos(/9' — 0) sinhr/T) /T'\ 



(A.ll) 



with C ■ 



(27r)3 



exp [/x/T]. Using the following identities for the modified bessel functions 



1 f'^'" 

In{z) = I-n{z) = — / exp [zees 0] cos(n0) d0 
27r Jo 

1 /"^ 

Kn{z) = K-n{z) = 7: exp [— zcosht] cosh(nt) dt 

2 J-oo 



and noting that 
f27r 



/ exp [z cos 9] sm{n9) d9 = 
Jo 



(A.12) 
(A.13) 

(A.14) 
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A Analytical derivations 



the rj- and ^-integrations can be done analytically: 

Vn =2Cpl j dp'T{p)p' j difjd^ 

r Qj- 
X cos(n0) Kif{ip)mT - ^o-^^Pt (/ cos{ip - (t>) + f sin(<y9 - 0)) 

X exp \pt cos((^ — 4>) sinh ryr] 
=2Cpo y dp'T{p)p' j dip j d^ (cos(n(P) cos(n(^) + sin(n<P) sin(n(^)) 

exp \pt cos(#) sinh ryr] 



dr 
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(A.15) 



In the first step, I used the identity of the K bessel function and substituted p = p' ■ pQ. 
In the second step, I substituted ^ = — cj) and expanded cos(n^) = cos(n((^ — #)). In 
the third step, I applied the trigonometric identities 



cos X cos?/ = - (cos(x — y) + cos{x + y)) 
sin X sin y = ^ (cos(x — y) — cos{x + y)) , 



(A.16) 
(A.17) 



A. 2 Integrals of the flow coefficients 
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so that I could solve the <P integral by using the identity for the / bessel function. In the 
last step, I used 



f'{<p) = f{^) 



sin(2(/?). 



(A.18) 



and again applied eq. (A. 17) to the sinus terms. Finally, the ^-integrated invariant yield 
at midrapidity reads 



AN 
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^exp [m/T] 47rp2 / d^p dp'T{p')p' 
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KiIof{ip)mT - KoIif{ip)pTdpT 
and the numerator of the elliptic flow is calculated as 

V2{Pt) = (^exp [n/T] AttpI ^ dipj^ dp'T{p)p' 
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dr 1 
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In the case of dpT = (radially constant freeze-out time), this simplifies to 
dN 
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C. 2 Deutsche Zusammenfassung der Diplomarbeit (German abstract to 
the thesis) 

Eine der grofien Fragen der Menschheit mit der sich Forscher der verschiedensten Diszi- 
plinen beschaftigt haben, ist die Frage nach den fundament alen Bausteinen der Natur bzw. 
ob es diese iiberhaupt gibt. Sie reicht zuriick bis in die Antike, zu dem Philosophen Empe- 
dokles, und den vier griechischen Elementen Erde, Wasser, Luft und Feuer. Im Mittelalter 
wurden immer mehr der heute bekannten chemischen Elemente klassifiziert und im Peri- 
odensystem gruppiert. Anfang des letzten Jahrhunderts stellte man fest, dass das Atom 
(griech. atomos: unteilbar) eine Substruktur, mit einem Kern aus Proton und Neutronen 
und einer Schale aus Elektronen, besitzt. Doch auch die Bausteine des Kerns offenbarten in 
inelastischen Streuexperimenten eine Substruktur. Die Fiille der neuen Teilchen genannte 
Hadronen, die bei diesen Experimenten erzeugt wurden, erforderte auch hier ein Ord- 
nungsschema. Aus dem "achtfachen Weg" von Murray Gell-Mann entwickelte sich dann 
das Quark-Modell, welches die Hadronen in Mesonen, bestehend aus einem Quark und 
einem Anti-Quark, und Baryonen, bestehend aus 3 Quarks, einteilt. 

Seit dem Beginn der Suche nach den elementaren Bausteinen der Natur benotigten die 
Forscher immer hohere Energien um die immer kleineren Langen aufzulosen. Nun suchen 
Physiker nach dem sogenannten Quark-Gluonen-Plasma (QGP), einem Zustand, in dem 
die Quarks nicht mehr eingesperrt sind in den Hadronen, sondern frei wie in einem Plasma. 
Mit der Kollision von schweren lonen am Relativistiv Heavy Ion Collider (RHIC), und bald 
mit dem Large Hadron Collider (LHC) am CERN, untersucht man den Ubergang zu dieser 
neuen Phase der stark wechselwirkenden Materie. 

Ein grofies Problem in diesem Zusammenhang ist die Frage, ob man iiberhaupt die 
Phasengrenze iiberschritten und ein QGP erzeugt hat, denn die Hadronisierung, also der 
Ubergang zuriick zur normalen Materie, findet auf extrem kurzen Zeitskalen statt. Daher 
konnen wir nur die Hadronen vom kinematischen Ausfrieren nachweisen, also die die keine 
Wechselwirkung mehr erfahren. Aber auch eine dynamische Beschreibung mithilfe der 
Quantenchromodynamik (QCD), der Eichtheorie der starken Wechselwirkung, ist nicht 
in naher Zukunft zu erwarten. Weder fiir die Kollision im Allgemeinen noch fiir den 
Vorgang der Hadronisierung im Speziellen. Somit sind wir auf phenomenologische Modelle 
angewiesen. 

Die Berechnungen innerhalb der pertubativen QCD (pQCD) mithilfe von gemessenen, 
parameterisierten Fragmentationsfunktionen sind sehr erfolgreich fiir Proton-Proton Kol- 
lisionen, aber in Schwerionenkollisionen mit einem dichten Phasenraum ist dieser Ansatz 

nur fiir hohe transversale Impulse anwendbar. Mit den ersten Resultaten von RHIC, 
die nicht mit pQCD zu erklaren waren, riickte ein phenomenologisches Modell namens 
Recombination in den Fokus. 

Im Gegensatz zu Fragmentation, bei der ein Quark aus einer harten Streuung mit hohem 

Impuls in viele Hadronen "zerfallt", die dann nur einen Bruchteil des Impulses tragen, ist 
die Idee bei der Recombination, dass in einem sehr dichten Medium ein Quark und ein 
Anti-Quark (drei Quarks) zu einem Meson (Baryon) recombinieren und sich dabei die 
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Impulse addieren. 

Eine der vielversprechensten Observablen bei der Suche nach dem QGP ist der so- 
genannte elliptische Fluss, welcher dem zweiten Fourierkoeffizienten V2 des im^arianten 
transversalen Impulsspectrums entspricht. Er quantifiziert die azimuthale Asymmetry des 
transversalen Impulses, welche aus der raumlichen Asymmetry in nicht-zentralen Kolli- 
sionen entsteht. Da der elliptische Fluss seinem Ursprung entgegenwirkt und damit ein 
selbstloschendes Verhalten zeigt, ist er damit hauptsaclilich von der friihen Pliase der 
Kollision gepragt und wiirde damit direkt vom Fluss der Quarks aus dem QGP abhangen. 

Der Erfolg der Recombination begriindet sich vorallem auf die Vorliersage einer uni- 
versellen Skalierung des elliptischen Flusses, wonach das Hadron V2 mit der Anzahl der 
Konstituentenquarks skaliert. Diese Vorhersage ist experimentell gut bestatigt und die 
Ergebnisse konnen die Daten von niedrigen bis mittleren transversalen Impulsen sehr exakt 
beschreiben. Wahrend das qualitative Verhalten des elliptischen Fluss dem Skalierungs- 
gesetz der Recombination folgt, hangen die quantitativen Ergebnisse von den zwei zusat- 
zlichen "Zutaten" ab: der Quarkdichteverteilung und der Hadronisierungshyperflache. 

In dieser Arbeit habe ich gezeigt, dass die einfliei^enden Modelle eine groi^e Wirkung in 
Bezug auf die Flusskoefhzienten haben. Fiir die Quarkverteilung benutze ich das Blast- 
Wave Modell, welches, basierend auf hydrodynamischen (Jberlegungen, ein transversal ex- 
pandierendes, lokal thermisches System beschreibt. Fiir die Parameterisierung des Druck- 
gradienten in niclit-zentralen Kollisionen geniigt ein simples geometrisches Argument um 
nicht nur den elliptischen Fluss V2, sondern audi den nachsten Koefhzienten V4 und deren 
relative Starke zu beschreiben. Ein Ergebnis des Blast-Wave Profils ist die beobachtete 
Massenskalierung von V2 und insbesondere ein negativer Fluss bei niedrigen transversalen 
Impulsen fiir schwere Hadronen und/oder hohe Schwerpunktsenergien. Seiche negativen 
Werte sind jetzt auch in den vorlaufigen J/i/j PHENIX-Daten zu sehen und scheinen 
damit die Anwendung der Blast-Wave Parameterisierung zu stiitzen. Das ist vorallem im 
Hinblick auf eine besonders iiberraschende Vorhersage in diesem Zusammenhang inter- 
essant. Denn Recombination in Verbindung mit dem Blast- Wave Modell prognostiziert 
eines abnehmenden (^2) fiir Schwerpunktsenergien oberhalb von ^/s = 5 TeV, welches 
moglicherweise am LHC beobachtet werden kann. 

Die Hadronisierung findet in einem 3-dimensionalen, zeitabhangigen Volume statt und 
fallt mit dem kinematischen Ausfieren zusammen, da die hadronische Streuphase nicht 
beriicksichtigt wird. Eine detaillierte Untersucliung der Hyperflache wiirde eine dynamis- 
che Behandlung erfordern, aber um grundlegende Einfliisse auf die Spektren zu unter- 
suchen ist dieser analytische Ansatz ausreichend, insbesondere well auch in der Hydrody- 
namik der Vorgang des Ausfrierens ein absolut nicht-triviales Problem ist. Wie meinen 
Ergebnissen zu entnehmen ist, haben eine raumliche Exzentrizitat beim Ausfrieren und 
die radiale Abhangigkeit der Ausfrierzeit den groi^ten Effekt auf die Flusskoeffizienten und 
besonders auf deren Verhaltnis ^4/(^2)^. Wahrend ein einfaches, zirkulares Ausfrieren 
dieses Verhaltnis um etwa einen Faktor 2 zu niedrig vorhersagt, erzeugt ein elliptisches 
Ausfrieren zusatzliche Beitrage zur Asymmetry beider Koefhzienten und kann damit die 
experiment ellen Daten sehr gut beschreiben. Und die radiale Abhangigkeit der Ausfrierzeit 
beeinflusst dann nochmal zusatzlich das Verhalten bei niedrigen transversalen Impulsen. 
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Die prasentierten Ergebnisse in meiner Arbeit sind ein klarer Hinweis darauf, dass 
Recombination der dominante Mechanismus zur Hadronisierung bei kleinen mittleren 
transversalen Impulsen ist, wobei die gute Beschreibung der Flusskoeffizienten und die 
Konstituentenquark-Skalierung ein starkes Zeichen fiir ein QGP in der friihen Phase der 
Kollision sind. Trotz der plienomenologischen Natur dieser analytisclien Untersuchung zur 
Hadronisierung in Schwerionenkollisionen durch Recombination, beschreiben die Ergeb- 
nisse nicht nur das qualitative Konstituentenquark-Skalierung der experimentellen Daten, 
sondern audi die quantitative Massenskalierung des elliptischen und hexadekupolen Flusses 
und das Verhaltnis der beiden Koeffizienten mit gro£em Detail. Aufgrund dieser guten 
tibereinstimmungen sollten die hervorstechenden Vorhersagen fiir den LHC detaillierter 
in dynamischen Untersuchungen betrachtet werden. 



